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Abstract 



In a previous paper Crater and Van Alstine applied the Two Body Dirac equations of constraint dynamics 
to quark-antiquark bound states using a relativistic extention of the Adler-Piran potential and compared 
their spectral results to those from other approaches which also considered meson spectroscopy as a whole 
and not in parts. In this paper we explore in more detail the differences and similarities in an important 
• subset of those approaches, the quasipotential approach. In the earlier paper, the transformation properties 

of the quark-antiquark potentials were limited to a scalar and an electromagnetic-like four vector, with the 
former accounting for the confining aspects of the overall potential, and the latter the short range portion. 
The static Adler-Piran potential was first given an invariant form and then apportioned between those two 
different types of potentials. Here we make a change in this apportionment that leads to a substantial 
improvement in the resultant spectroscopy by including a time-like confining vector potential over and above 
^ l' the scalar confining one and the electromagnetic-like vector potential. Our fit includes 19 more mesons than 

O , the earlier results and we modify the scalar portion of the potential in such a way that allows this formalism 

^ • to account for the isoscalar mesons rj and 77' not included in the previous work. Continuing the comparisons 

of formalisms and spectral results made in the previous paper with other approaches to meson spectroscopy 
, we examine in this paper the quasipotential approach of Ebert, Faustov, and Calkin. 

> 

o 

00 

^ ! 1 Introduction 
^' 

, There are a number of strategies in computational treatments of quantum chrome dynamics that emerge in the 

^■"■^ ' study of meson spectroscopy. One is to set up a discrete lattice analog of the full quantum field theory. A 

second is to first make analytic approximations which replace the quantum field theoretic problem by a classical 
variational problem involving an effective Lagrange function and action. The latter approach has been exploited 
^ , by Adler and Piran [J and in a previous paper Crater and Van Alstine gave a detailed account of applications 

^ ' of the Two Body Dirac equations (TBDE) of constraint dynamics to the meson quark-antiquark bound states 

[2] using a relativistic extension of the Adler-Piran potential. That paper also included a comparison of this 
approach to others [3]- [6] who, like ours, considered the whole spectrum instead of just selected parts. 

Here we update the results presented in [2] in four ways. First, we include 19 more mesons not included 
in the previous work. Second, we obtain a substantial improvement in our fit to most all of the mesons by 
allowing the confining interaction, pure scalar in [5], to take on a timelike vector portion. We still include the 
electromagnetic-like vector potential used previously. Third, we extend the relativistic Schrodinger-like form of 
the TBDE to include isoscalar mixing, thus incorporating the isoscalar mesons rj and 77'. And finally we critically 
examine, by comparison with the TBDE, aspects of quasipotential approaches including a recent one presented 
in [7] as well as in [5]. 

In Sec. 2 we give a short review of the relativistic two-body constraint formalism, distinguishing between our 
new approach used for confining given in this paper and the one presented in [5] and including a discussion of 
the closely related quasipotential approach. In Sec. 3 we review the static Adler-Piran potential and how we 
apportion it between the three invariant potential functions A(r), S{r), and V{r) used in our TBDE. In Sec. 
4 we present our main new results on meson spectroscopy. In Sec. 5 we include our treatment of SU(3) mixing 
and in Sec. 6 we discuss the meson spectral results of [7J including the advantages and shortcomings of their 
quasipotentials bound state formalism. 
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2 Review of Relativistic Two-Body Formalisms 



2.1 Two-Body Constraint Approach 

When the interaction and the masses are known, a common starting point in describing the relativistic two-body 
bound state problem is the Bethe-Salpeter equation [8]. The Bethe-Salpeter equation is, however, usually not 
considered in its full four-dimensional form due to the difficulty of treating the relative time coordinate [9]. 
Numerous truncations of the Bethe-Salpeter equation have been proposed for the relativistic two-body problem 
[TUl [TT] . Some of these types of approximate methods have previously been applied with considerable success 
to the qq meson spectrum |12|-|18), [2]. The ladder approximation and the instantaneous approximation of the 
Salpeter equation have been widely used. It should be noted, however, that the simple ladder approximation 
and the Salpeter equation do not lead to the correct one-body limit [12], and do not respect gauge invariance 
[20) . Crossed-ladder diagrams must be included to insure gauge- invariant scattering amplitudes. 

The Two Body Dirac equations of constraint dynamics provide a covariant three-dimensional truncation of 
the Bethe-Salpeter equation. Sazdjian [21] has shown that the Bethe-Salpeter equation can be algebraically 
transformed into two independent equations. The first yields a covariant three-dimensional eigenvalue equation 
which for spinless particles takes the form 

nio + n2o + 2<i>^}j^(xi,x2) = o, (1) 

where Hio — pf + mf . The quasipotential is a modified geometric series in the Bethe-Salpeter kernel K such 
that in lowest order in K 

= TriwS{P ■ p)K, (2) 

where P = pi + P2 is the total momentum, p — 772P1 ~ ''I1P2 is the relative momentum, w is the invariant total 
center of momentum (cm.) energy with — —tip'. The rj^ must be chosen so that the relative coordinate 
X = xi — X2 and p are canonically conjugate, i.e. ?7i + 7?2 = 1- The second equation overcomes the difficulty of 
treating the relative time in the center of momentum system by setting an invariant condition on the relative 
momentum p, 

(Hio-H2o)*(a;i,X2) = = 2P-p*(xi,a;2). (3) 

Note that this implies p'^* = p^^^ = {rif^" + P'^P'')p^'^ in which = Pf^/w is a time like unit vector (P^ = -1) 
in the direction of the total momentum. 

One can further combine the sum and the difference of Eqs. ([T]) and ^ to obtain a set of two relativistic 
equations one for each particle with each equation specifying two generalized mass-shell constraints 

Hi'^ixi,X2) = [pf + + 'Pni)'^{xi,X2) = 0, i = 1,2, (4) 

including the interaction with the other particle. These constraint equations are just those of Dirac's Hamiltonian 
constraint dynamic^ [551 [53]. In order for the two simultaneous wave equations of Q to have solutions other than 
zero, Dirac's constraint dynamics stipulate that these two constraints must be compatible among themselves, 
['Hi,'H2]^' = 0, that is, they must be first class. With no external potentials, the coordinate dependence of the 
quasipotential would be through x and the compatibility condition becomes [pf —P2, $u)]^ = P^d^^/dx'^ = 
0. In order for this to be true in general, must depend on the relative coordinate x only through its component, 
x_L, perpendicular to P, 

a;^[ = (,7''''-fP'^P'')(xi-X2).. (5) 

Since the total momentum is conserved, the single component wave function in coordinate space is a product 
of a plane wave eigenstate of P and an internal part tp [24j , depending on this x± H 

We find a plausible structure for the quasipotential $i„ by observing that the one-body Klein-Gordon equation 
(p^ + ■m'^)ip — (p^ — + nP)ip = takes the form (p^ ~ e"^ + m? + 2mS + + 2eA — A'^)'ip = when one 
introduces a scalar interaction and timelikc vector interaction via m ^ m + S and e ^ e — A. In the two- 
body case, separate classical ^ and quantum field theory [53] arguments show that when one includes world 

1 These equations were originally proposed in the form of classical generalized mass shell first class constraints "Hi = (p^ + + 
$i) ~ 0, and their quantization 'Hi'i' = without reference to a quantum field theory. For the classical Hi to be compatible, their 
Poisson bracket with one another must either vanish strongly or depend on the constraints themselves, {^1,^2} ~ 0. The simplest 
solution of this equation is = $2, a kind of relativistic third law condition, together with their common transverse coordinate 
dependence ^w{x±), just as with its quantum version. 

^We use the same symbol P for the eigenvalue so that the w dependence in Eq. ^ is regarded as an eigenvalue dependence. The 
wave function can be viewed either as a relativistic 2-body wave function (similar in interpretation to the Dirac wave function) or, 
if a close connection to field theory is required, related directly to the Bethe Salpeter wave function x by |21| 'J = —m5{P ■p)'HwX = 
-■KiS{P ■ p)'H20X- 
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scalar and vector interactions then depends on two underlying invariant functions S{r) and A{r) through 
the two-body Klein-Gordon-like potential form with the same general structure, that is 

'i>ni ^2m,,S + + 2eniA- A^. (6) 

Those field theories further yield the cm. energy dependent forms 

rrini = mim2/w, (7) 

and 

Etu = (w^ ~ ml ~ ■ml)/2w, (8) 

ones that Tododov [33] introduced as the relativistic reduced mass and effective particle energy for the two-body 
meson system. Similar to what happens in the nonrelativistic two-body problem, in the relativistic case we have 
the motion of this effective particle taking place as if it were in an external field (here generated by S and A). 
The two kinematical variables ([7]) and ^ are related to one another by the Einstein condition 

sl-ml = b\w), (9) 

where the invariant 

b^iw) = {w^ - 2w'^{ml + m\) + (m\ - to^)^)/4u;^ (10) 
is the cm. value of the square of the relative momentum expressed as a function of w. One also has 

}?{w)=E\-m\^e\-m\, (11) 

in which e\ and £2 are the invariant cm. energies of the individual particles satisfying 

El + £2 = S\ — ei = {ml — m\)/w. (12) 

In terms of these invariants, the relative momentum appearing in Eq. ^ and Q is given by 

p^^ = {e2Pt - ew^^)lw, (13) 

so that rji+rj2 = (ei +£2)/^ = 1- In [22] the forms for these two-body and effective particle variables are given 
sound justifications based solely on relativistic kinematics, supplementing the dynamical arguments of |25l and 

m- 

Originally, the Two Body Dirac equations of constraint dynamics arose from a supersymmetric treatment 
of two pseudoclassical constraints (with Grassmann variables in place of gamma matrices) which were then 
quantized |28j-|31j. Sazdjian later derived j2l) different forms of these same equations, just as with their spinless 
counterparts above, as covariant and three-dimensional truncation of the Bethe-Salpeter equation. The forms 
of the equations are varied (see Appendix A), but the one that is the most familiar is the "external potential" 
form similar in structure to the ordinary Dirac equation. For two particles interacting through world scalar and 
vector interactions they are 

Slip EE 75i(7i • (pi - Ai) + mi + S-i)* = 0, 

S21P = 752(72 • iP2 - A2) + m2 + 52)* = 0. (14) 

Here is a 16 component wave function consisting of an external plane wave part that is an eigenstate of P 
and an internal part ip = ipix^). The vector potential A'^ was taken to be an electromagnetic-like four- vector 
potential with the time and spacelike portions both arising from a single invariant function A. The tilde on 
these four- vector potentials as well as on the scalar ones Si indicate that they are not only position dependent but 
also spin-dependent by way of the gamma matrices. In this paper we allow for the presence of a timelike portion 
arising from an independent invariant function V. In either case, the operators Si and ^2 must commute or at 
the very least [Si,S2]ip = since they operate on the same wave function. This compatibility condition gives 
restrictions on the spin dependence which the vector and scalar potentials 

if = if (A(r), V{r),p^,P, w, 7i, 72), = S,{S{r), A{r),p^,P, w, 7i, 72)- (15) 



^In a perturbative context, i.e. for weak potentials, that would mean that this aspects of Af is regarded as arising from a 
Feynman gauge vertex coupling of a form proportional to 7^72^^ (^^e Appendix A). 

*In a perturbative or weak potential context, that would mean that this aspect of Af is regarded as arising from an additional 
vertex coupling proportional to —-fi ■ P ^2' Similarly in a perturbative or weak potential context Si is regarded as arising from 
a vertex coupling proportional to I1I2S'. (See Appendix A). 

^The 75 matrices for each of the two particles are designated by 75^ i = 1,2. The reason for putting these matrices out front of 
the whole expression is that including them facilitates the proof of the compatibility condition, see |28| , I24| . 
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are allowed to have[f| in addition to requiring that they depend on the invariant separation r = \px^ through the 
invariants A(r), V(t)^ and S(t) . The covariant constraint ([3]) can also be shown to follow from Eq. (|14p. We 
give the explicit connections between A^, Si and the invariants A(t)^ ^ {t\ Sir) in Appendix A. The Pauli 
reduction of these coupled Dirac equations lead to a covariant Schrodinger-like equation for the relative motion 
with an explicit spin-dependent potential 

[v\ + ^^{A(r),V{T),S{T),p^,P,w,au<y2)^i'^^h'{w)^^ , (16) 

with }?{w) playing the role of the eigenvalue^ This eigenvalue equation can then be solved for the four-component 
effective particle spinor wave function -0^ related to the 16 component spinor ^'(^^-l) in Appendix A. 

The set of equations ()14p and the equivalent Schrodinger-like equation (|16p possess a number of important 
and desirable features. First, they reduce to the correct one-body Dirac form when one of the two constituents 
becomes very massive. (The Salpeter equation does not have this important property.) Second, the generalized 
three-dimensional Schrodinger equation ([T5|) is quite similar to the nonrelativistic Schrodinger equation and it 
indeed goes over to the correct nonrelativistic Schrodinger equation in the limit of weak binding. One can thus 
employ familiar techniques to obtain its solutions. Third, Eq. (jl6p can be solved nonperturbatively for both 
QED bound states (e.g. positronium and muonium) and QCD bound states (i.e. bound states obtained from 
two-body relativistic potential models for mesons) since every term in is nonsingular in the sense that they 
are less attractive than — l/4r^ (no delta functions or attractive l/r^ potentials for example). Thus, unlike with 
the 1/m non-relativistic and semirelativistic expansions, the covariant Dirac formalism itself introduces natural 
cutoff factors that smooth out singular spin-dependent interactions, there being no need to introduce them by 
hand (- see [2H] , [22] , [29] , [33] and Sec. 6.4) as in other approaches. Fourth, the relativistic potentials appearing 
in these equations are directly related through Eq. ([2]) to the interactions of perturbative quantum field theory, 
while for QCD bound states they may be introduced semiphenomenolgically through A(r) and Sir) (and in this 
paper V(r)). Fifth, these equations have been tested analytically [34, and numerically [22] against the known 
perturbative fine and hypcrfinc structures of QED bound states and related field theoretic bound states. The 
(nonperturbative) successes with these QED bound states provide strong motivation for applying the constrained 
Dirac formalism to meson bound states as in [2 . Sixth, these equations provide a covariant three-dimensional 
framework in which the local potential approximation consistently fulfills the requirements of gauge invariance 
in QED |35| . Finally, the same general structures of the Darwin, spin-orbit, spin-spin and tensor terms in 
<f>tt,(A, y, 5) of Eq. p6]l and p5)) . responsible for the accurate hyperfine structures of QED bound states arising 
naturally from the TBDE formalism when A = —a/r and S — V — are used with the only alteration in its 
application to the QCD bound states being that A^ V, S are apportioned appropriately from the Adler-Piran 
potential. 

We emphasize the importance of the nonperturbative QED bound state test. (See Appendix C for a review 
of the application of the TBDE to QED.) Many of the wave equations used in the standard approaches to 
QED bound states have been modified to include QCD inspired potentials and then applied to nonperturbative 
numerical calculations of QCD bound states without first testing those approaches nonperturbatively in QED. By 
this we mean that the accepted perturbative results of those equations (QED spectral results correct through order 
a*) have not been replicated using numerical methods. Somnierer et al. [4] have shown that the Blankenbecler- 
Sugar equation and the Gross equations fail this test. This indicates a danger in applying such three-dimensional 
truncations of the Bethe-Salpeter equation: if failure occurs in their applications to QED bound states this 
brings into question the spectral results of similar nonperturbative (i.e. numerical) approaches based on the 
same truncations when applied to QCD bound states. This would be true especially when the only difference 
between the vector portions of the QED and QCD potentials would be the replacing of the QED —a/r by a 
similar A{r) from QCD. 

In ^ we presented details of the application of this formalism to meson spectroscopy using a covariant 
version of the Adler-Piran static quark potential. Note especially that the equations used there displayed a 
single <i>(yl(r), S{r),p±, P, w, ci, tT2, ) in Eq. (fTH)) . It depends on the quark masses through factors such as those 
that appear in Eq. ([6]). However its dependence is the same for all quark mass ratios - hence a single structure 
for all the QQ, qQ, and qq mesons in a single overall fit. We found that the fit provided by the TBDE for the 
entire meson spectrum (from the pion to the excited bottomonium states) competes with the best fits to partial 
spectra provided by other approaches and does so with the smallest number of interaction functions (just A{r) 

^The dependence of the scalar potentials Si on the invariant A(r) responsible for the electromagnetic-like potential is seen in I24| 
and I2bl to result from the way the scalar and vector fields combine. That combination without the presence of the independent 
time- like portion leads to a two-body Klein-Gordon-like potential portion of <I>u, to be of the form given in Eq. (O. 

^In the presence of an additional and independent time-like vector interaction V, we assume the scalar and vector fields combine 
in such a way that leads to a two-body Klein-Gordon- like potential portion of of the form 2mu, 5-1-5^-1- 2e-w A — -\- 2e-w V — 
instead of that given in Eq. {G} . 

*Due to the dependence of on ui, this is a nonlinear eigenvalue equation. 
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and S{r)) without additional cutoff parameters necessary to make those approaches numerically tractable. We 
also found that the pion bound state displays some characteristics of a Goldstone boson. That is, as the quark 
mass tends to zero, the pion mass (unlike the p and the excited tt) vanishes, in contrast to almost every other 
relativistic potential model. (For more discussion on this see footnote 23 below). 



2.2 Two Body Quasipotential Approaches 

Also presented in [2] was a detailed comparison between the meson spectroscopy results of our model and 
those of several other approaches: one based on the Breit equation [S], two on truncated versions of the Bethe- 
Salpeter equation [3] , [1] , and one on a quasipotential approach [S] . We explore in this section in more detail the 
differences and similarities between our approach and the quasipotential approach. The quasipotential equation 
was first introduced by Logunov and Tavkhelidze [TO]. In its homogeneous form, that equation describes a 
two-particle relativistic composite system with its cm. momentum space form (in the notation used here p is 
the relative momentum given in Eq. ()13|) ) for spinless particles given by 



(u; - ^p2 + mf - ^p2 + m2)^'^(p) =. J \/(p, q,zi;)^f^(q) (17) 

where 5'ii,(p) is the quasipotential wave function projected onto positive- frequency states and V{p,q,w) is the 
quasipotential calculated by means of the off-energy shell scattering amplitude (so that the respective cm. 
energies of the two particles are not given by the above square roots but by Eq. (|12|) ). The corresponding 
inhomogeneous quasipotential equation is of the general form 

T{p,ci,w) + Vip,q,w)+ J V{p,k,w)G^{k)V{k,ci,w) ^0, (18) 

a linear integral equation of the Lippmann-Schwinger type relating the quasipotential to the off-energy shell 
extrapolation of the Feynman scattering amplitude. The choice of this equation and the accompanying ho- 
mogeneous equation is not unique |11| . For example, the Green function Giu(k) has only its imaginary part 
determined by requiring the condition of elastic unitarity on Eq. ([18]). Todorov 10 took advantage of this 
nonuniqueness to write down a local version of the corresponding homogeneous equation of the form 

(p2 - b') 0(P)+^ / K„(p,k)0(k)^ = 0, (19) 

with (/)(p) the wave function in momentum space. In |32j Crater and Van Alstine showed that the spinless version 
of Eq. (|16|) in the case of QED {V = S = 0) has, in the cm. frame, the form (see Eq. (66b) and discussion 
below Eq. (73b) of that paper) 

(p2 _ _ ^)2 + + 1 ^ 1 (vf;)2)^ = 0, (20) 

where 

g = InG, 

^ = (l-2A/u;)i/2- (21) 

As was pointed out in that paper, for A — —a/r, this reduces for weak potentials to the minimal or gauge 
structure form postulated by Todorov, 

[(p-A)'-(£„-A")2+m^]V = 0. (22) 

Although not noticed at the time, Eq. (|20]) does in fact additionally have this minimal structure not only for 
arbitrary strength couplings but also for potentials not restricted to Coulomb potentials, provided just that 

A" = A, 

A = -^V5/„ (23) 
where Ig is the space reflection operator satisfying 

IJ{r) =/(-r). (24) 



'For hermitian potentials, that condition has the symbolic form of T — Tt = T^(G — G^)T. 
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Aneva, Karchev, and Rizov [35] developed the weak potential version of this two-body Klein-Gordon equation 
for two spin combinations: for one spin-zero and one spin-one half particle and for a spin-one-half particle- 
antiparticle pair. For the latter it has the form, 

(p2 - b') ^,^,^(p)+^ I uy^{p)v,,i-k)Y^ip,k)vy^{-p)u,Amy,y,{^)-^ = 0- (25) 

Expressing the on shell free four-component Dirac spinors in terms of two-component Pauli spinors and assuming 
the local quasipotentials 

V„(p, k) =V„(p - k) =^(p - k)75'72^ + V(p - k)/3i/32 + 5(p - k)lil2, (26) 

Eq. (j25p can be brought to a four-component wave equation form superficially similar to Eq. (jl6p in the cm. 
frame. We write it as 

[v" + V^(A(r), ^(r), 5(r), p, <ti, ^3)^ 0(r) = \?(w)ci,{v). (27) 

However, there are distinct differences. First of all, the spin structure of $ui is not identical to that of even 
if the functions A(r)^ ^(^)i S(r) are the same. The reason is that the spin dependence of the vector and scalar 
potentials and Si and in particular the minimal type of context in which they appear in Eqs. (|14p arise 
from (nonlinear) hyperbolic functions (see Appendix A and [37|, [SQ) of matrices such as appear in Eq. ([26]) . 
Without that hyperbolic structure the external potential forms in which the minimal structures appear would 
be absent. To reproduce the effects of those nonlinear functions in Eq. (|27p one would have to supplement 
Eq. (|26|) with types of invariants other than just scalar and vector [39], [40], [26] (a pseudovector invariant for 
example). Secondly, the desirable minimal scalar structures as appear in the first line on the right hand side 
of Eq. (1351) below would not appear in V^, without including higher order diagrams that would again require 
invariants other than just scalar and vector. These minimal scalar structures are not only desirable, they arise 
naturally and strictly from 0(l/c^) expansions of classical and quantum field theoretic potentials [25, [26] and 
from gauge invariance considerations (see Todorov in [10], and also [36], [41], and p4|). In a later section we 
discuss the recent work of [7], which uses a quasipotential equation similar to Eq. p7p in meson spectroscopy 
calculations. 



3 The Adler-Piran Potential for the Two Body Dirac Equations. 

In this section we use the relativistic Schrodinger-like Eg. dTB)) to construct a relativistic naive quark model by 
choosing the three invariant functions A, V and S to incorporate the Adler-Piran static quark potential [T] . 
This potential was originally obtained from the QCD field theory through a nonlinear effective action model 
for heavy quark statics. Adler and Piran used the renormalization group approximation to obtain both total 
fiux confinement and a linear static potential at large distances. Their model uses nonlinear electrostatics with 
displacement and electric fields related through a nonlinear constitutive equation with the effective dielectric 
constant given by a leading log-log model which fixes all parameters in their model apart from a mass scale A. 
Their static potential also contains an unknown "integration constant" \Jq in the final form of their potential 
(hereafter called VAp(f')). We insert into Ea.([T6| invariants A^V and S with forms determined so that the 
sum A + V ^ S appearing as the potential in the nonrelativistic limit of our equations becomes the Adler-Piran 
nonrelativistic QQ potential (which depends on two parameters A and fJo) plus the Coulomb interaction between 
the quark and antiquark. That is, 

VAp(r) + Vcoui = A(;7(Ar) -^U^) ^ — = A + V + S . (28) 

r 

As determined by Adler and Piran, the short and long distance behaviors of U (Ar) generate known lattice and 
continuum results through the explicit appearance of an effective running coupling constant in coordinate space. 
That is, the Adler-Piran potential incorporates asymptotic freedom through 

At/(Ar << 1) - l/(rlnAr), (29) 

and linear confinement through 

A[/(Ar >> 1) - AV. (30) 

In addition to obtaining these leading behavior analytic forms for short and long distances, they converted the 
numerically obtained values of the potential at all distances (short, intermediate, and long distances) to compact 
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analytic expressions. The explicit closed form expressions [T* for U{Ar) are different for each of the four regions, 
and are linked continuously. Letting x = Ar, 



U{x) = -(167r/27)(l + aia;"2)/(wp)/wp, < a: < 0.0125, 

Wp = l/iasx)'^, 

U{x) = X + a(a;/0.125)^, 0.0125 < x < 0.125, 

E = f3 + ^\n{l/x)+5{\n{l/x)f +e{\n{l.x)f 

U{x) = K' + a'\TLX + l3'{\TLxf +-f'{lnxf + 5'{\nxf + e'{lnxf, 0.125 < Ar < 2, 

U{x) = k[cix + C2\nx + + — + c^), 2 < Ar < CX3. (31) 

\/x X 

The function / is defined by = /(ln/ + Clnln/); C = 2(51 - 19n//3)/(ll - 2n//3)2 = 64/81 for n/ = 3 . The 
various constants ai to as, if, a, /3, 7, 5, e, if', a', 7', (5', e',and ci to C5 are given by the Adler-Piran leading 
log-log model [1] and are not adjustable parameters. We modify these closed forms so that the connections 
between different regions are continuous in second derivatives. The nonrelativistic analysis used by Adler and 
Piran, however, does not determine the relativistic transformation properties of the potential. How this potential 
is apportioned between vector and scalar is therefore somewhat, although not completely, arbitrary. 

In earlier work [13] we divided the potential in the following way among three relativistic invariants A, V 
and 5' for all x = Ar. (In our former construction, the additional invariant V was responsible for a possible 
independent timelike vector interaction.) 

S = 77(A(cix + C2 ln(.T) + C5 -f ?7o), 

V X 

V ={l- f])A{cix + C2 ln(x) + + C5 + Uo), 

^x 

A = L/(x)-A(cix-|-C2ln(x)-|--^-f eg), (32) 

^x 

in which 77 = ^. That is, we assumed that (with the exception of the Coulomb-like term (€4/0;)) the long distance 
part was equally divided between scalar and a proposed timelike vector. 

In the present investigation, we compute the best fit meson spectrum for the following apportionment of the 
Adler-Piran potential: 

A = exp(-/3Ar)[y^p-^] + ^ + ^, 

r r r 

V + S = yAP + — -A=CK4p--)(l-exp(-/3Ar))=i/, (33) 
r r 

In order to covariantly incorporate the Adler-Piran potential into our equations, we treat the short distance 
portion as purely electromagnetic-like (in the sense of the transformation properties of the potential). The 
attractive (04 = —0.58) QCD-Coulomb-like portion (not to be confused with the electrostatic Vcoui = eie2/r) is 
assigned completely to the electromagnetic-like part A. That is, the constant portion of the running coupling 
constant corresponding to the exchange diagram is expected to be electromagnetic-like (~ 71^72)- Through 
the additional parameter /3, the exponential factor gradually turns off the electromagnetic-like contribution (i.e. 
A) to the potential at long distance except for the 1/r portion mentioned above, while the scalar and timelike 
portions (i.e. S and gradually turn on, becoming fully responsible for the linear confining and subdominant 
terms at long distance. We choose not to consider an apportionment function with a large number of parameters 
as the simple exponential gives a single length scale for the turning of the potential from electromagnetic-like to 
scalar and timelike. Altogether our three invariant potential functions depend on three parameters: A, J/q, and 
/?. We furthermore let a free parameter ^ divide the relative portions of U as follows 

S = eZ^-e(^AP"-)(l-exp(-/3Ar)), 
r 

V = U-S = {l-i){VAp--){l-eM~P^r)). (34) 

r 

This differs from the division in the earlier work 13:. Also, the earlier work did not include the effects of the 
tensor interaction or spin-orbit difference terms or the u — d quark mass differenced^ (see Eq. (I55t below). In 



'^'^In the present treatment, we treat the entire interaction present in our equations, thereby keeping each of these effects, fn our 
former treatment |13| we also performed a decoupling between the upper-upper and lower-lower components of the wave functions 
for spin-triplet states which turned out to be defective but which we subsequently corrected in our numerical test of our formalism 
for QED [^. 



7 



[2] , Crater and Van Alstine chose ^ = 1 and thus assumed that the scalar interaction is solely responsible for 
the long distance confining terms. 

When inserted into the constraint equations, V, S and A become relativistic invariant functions of the invari- 
ant separation r = ^/x^ . The covariant structures of the constraint formalism then automatically determine the 
exact forms by which the central static potential is supplemented with accompanying relativistic spin-dependent 
and recoil terms. 

4 Meson Spectroscopy from the Schrodinger-like form of the Two- 
Body Dirac Equations 

4.1 Center of Momentum form of Covariant Pauli-Schrodinger Reduction of the 
Two-Body Dirac Equations 

In Appendix A we outline the steps needed to obtain the explicit cm. form of Eq. (fT6|). That form is j42j . 

001, m, 

{p^ + $(r,mi, m2, 0-1,0-2)}^+ = 
={p2 + 2m^S + S^ + 2e^A - + 2e^V -V^ + <^d 

+ L • {(T^+(T2)'^SO + (Ti-fcr2-fL • {(T^+(T2)^SOT 
+ Cri Cr2<l>SS -I- (3(Ti fcr2 f - (Ti (T2)<f>T 

+ L • (cr^-<T2)$soD + «L • (TiX(T2$sox}^/'+ 

= b^^^+. (35) 

The detailed forms of the separate quasipotentials $i are given in Appendix A together with their forms for 
weak potential and in the static limit. In Appendix B we give the radial forms of Eq. ([55)1 . The subscripts of 
most of the quasipotentials are self explanatory. After the eigenvalue 6^ of ([BS)) is obtained, the invariant 
mass of the composite two-body system w can then be obtained by inverting Eq. (|10|) . It is given explicitly by 

w = \Jh'^ + mj + \Jb'^ + ml. (36) 

The structure of the linear and quadratic terms in Eq. psp as well as the Darwin and spin-orbit terms, are 
plausible in light of the discussion given above Eq. ([6]), and in light of the static limit Dirac structures that 
come about from the Pauli reduction of the Dirac equation (see Eqs. (ffS]) below for the two-component Pauli 
reduction of the Dirac equation). Their appearances as well as that of the remaining spin structures are direct 
outcomes of the Pauli reductions of the simultaneous TBDE Eq. ([M]). 

4.2 Spectral Results 

Theoryl (abbreviated by Thl) has two invariant interaction functions: A{r) for the short distance behavior and 
S{r) for scalar confinement. Theory 2 (abbreviated Th2) has three invariant interaction functions including 
the previous two plus V{r) for timelike vector confinement. In Appendix A.l we outline how these invariant 
interaction functions are related to what we call vertex invariants {J^{r), G{r), C{r)), which define the covariant 
gamma matrix interaction structures which enter into the hyperbolic form of the TBDE as seen in Eqs. (jA.2|lAT3| 
IA.7IIA.12]) . and the related energy and mass potentials Ei^2 and Afi,2, which characterize the external potential 
forms of the TBDE as seen in Eqs. (|A.13IIA.16p . The distinction between Thl and Th2 includes more than the 
partial cancellation — = iY^(2^ — 1) between the quadratic scalar and timelike vector interactions. It 
also includes their partial cancellations for spin-orbit and Darwin terms (see Appendix A. 4). In this section we 
present spectral results with (Th2) and without (Thl) the added timelike invariant function V{r). 

We display our results in Tables [TIIT51 The first table lists the best fit values for the quark masses and the 
potential parameters A, AUq, (3 for Thl (scalar only confinement) and Th2. The ratio that optimized the fit 
for Th2 is ^ = 0.704, (see Eq. (p4|) '). In the first two columns of Tables [2][9l we list quantum numbers and 
experimental rest mass values (in GeV) and experimental errors listed parenthetically (in MeV) for 105 known 
mesons. We include all well known and plausible candidates listed in the standard reference ([44]). We omit 
only those mesons with substantial flavor mixing, like the and 77' mesons . In the tables, the quantum numbers 

^^The subscript on quasipotential #13 refers to Darwin. It consist of what are called Darwin terms, those that are the two-body 
analogue of terms that accompany the spin-orbit term in the one-body Pauli reduction of the ordinary one-body Dirac equation, 
and ones related by canonical transformations to Darwin interactions 1251 1431 . momentum dependent terms arising from retardation 
effects. The subscripts on the other quasipotentials refer respectively to SO (spin-orbit), SOD (spin-orbit difference), SOX (spin- 
orbit cross terms), SS (spin-spin), T (tensor), SOT (spin-orbit-tensor) 
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PcirciiiiGtcr 
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TVi9 
± IIZ 


mi, 




A OKI CoAf 


rric 


i.04D KjG V 


1 PoA/' 

i.ooo V 


rris 


n 9S7/I PoA/' 
U.Zo < 4 KjG V 


u.ou / y KjG V 


rriu 


n 071 r'o\f 
(j.v 1 lo V 


u.uyoo KjG V 


fnd 


n n77i r^^iA/' 

U.U ( ( i KjG V 


u. iU4t) (^e V 


A 


n 991 CnV 


n 99^^ Pr.A/^ 


A TT 
AUq 


i.oio V 


1 770 PnAA 


6 


1.502 


4.408 




1 


0.704 



Table 1: Parameters for Theory 1 and 2. The 4 potential parameters A, AUo, l/(/3A), and ^ are respectively the 
QCD scale factor, the Adler-Piran integration constant, the vector-scalar transition distance, and the S/ {S + V) 
ratio. 

listed are those of the ip^ part of the 16-component wave function. In the third and fourth columns are the 
theoretical values, with Thl referring to the results without the timelike vector interaction and Th2 with the 
timelike vector interaction. In the fifth and sixth columns we give the differences between our theoretical results 
and the experimental and in the last two columns the contributions of each theoretically computed values to the 
total of 237 for Thl and 173 for Th2El- 

To generate the fits, in addition to varying the five quark masses we vary the parameters A, Uq, /3, and ^ in 
the apportioned static Adler-Piran potential in A, V, and S. Those invariants are put into our relativistic wave 
equations just as we have inserted the invariant Coulomb potential A = —a/r (but with = S* = 0) to obtain the 
results of QED bound states [Ml HI]. Note especially that we use a single S) for Thl and a single V, S) 
for Th2 for all quark mass ratios. Hence in each theory we use a single structure for all the QQ, qQ, and qq 
mesons in a single overall fit. The entire confining part of the potential transforms as a world scalar for Thl and 
combined timelike and scalar for Th2. Since ^ > 0.5 in our equations, this structure leads in both models to 
linear confinement at long distances and quadratic confinement at extremely long distances (where the quadratic 
contribution outweighs the linear term 2mwS in Thl and S'^ — V^ outweighing the linear terms 2myjS -\-2ewV 
in Th2). At distances at which exp(— /3Ar) << 1, the corresponding fine and hyperfine structures producing 
spin-orbit, Thomas, Darwin, spin-spin, and tensor terms (the last two are relatively small in that domain) 
are dominated by the confining interaction, while at short distances (exp(— /3Ar) ~ 1) the electromagnetic-like 
portion of the interaction gives the dominant contribution to the fine and hyperfine structures. Furthermore 
because the signs of each of the spin-orbit and Darwin terms in the Pauli-form of our TBDE are opposite for the 
scalar and vector interactions (see Appendix A. 4) , the spin-orbit contributions of those parts of the interaction 
produce opposite effects with degrees of cancellation depending on the size of the quarkonia atom. Another 
point to make is that because of the various sizes of the quarkonia atoms and the cm. energy dependence the 
behavior of is sharply different for the light mesons compared with the heavy ones. This may possibly 
account for the ability of our formalism to obtain good fits for the light meson hyperfine splittings while at the 
same time giving good overall fits to the heavy mesons. 

We obtain the meson masses given in columns three and four as the result of a least squares fit using the 
known experimental errors from the Particle Date Group (PDG) tables j44j and an assumed calculational error 
of 1.0 MeV. We employ the calculational error not to represent the uncertainty of our algorithm but more to 
prevent the mesons that are stable with respect to the strong interaction from being weighted too heavily. Our 

is per datum (105) minus parameters (8 or 9). In Table [TJ the value of /? for Thl implies that (in the best 
fit) as the quark separation increases, our apportioned Adler-Piran potential switches from primarily vector to 
scalar at about (/3A)~^ ~0.60 fermi. This shift is a relativistic effect since the effective nonrelativistic limit 
of the potential {A + S) exhibits no such shift (i.e., by construction /? drops out). For Th2, this distance is 
substantially less, (^A)"! -0.22 fe rmi. 



Table [To] lists the 8 meson families, their respective x^ contributions and their averages. The most striking 
feature is that as the quark masses increase from the lightest to the heaviest, the differences of the respective x^ 

^■^The reader of ^ may notice that the total of the model there of 101 (corresponding to Thl here) was substantially lower 
than the 237 that we found here. There are several reasons for this difference. The main reason is that in this paper we do not 
include a 5% addition to the calculational uncertainty based on the total meson widths. Also, there are 19 more mesons in the 
present model, many of which were difficult to fit. Thirdly, the experimental errors changed. Fourth, many of the newer mesons 
(for example the rif,) not only added more to the x'^ from their own fit, but also indirectly to the older ones (e.g. the 1^S\ T meson). 
There is no difference in the parameter sets and potentials in Thl and those used in [2) although the values are different. 
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ud 


mesons 


Exp. 


Thl. 


Th2. 


Exp.-Thl. 


Exp.-Th2. 


X"-Thl. 


X'-Th2. 


TT : 


ud 1 ^So 


0.140(0.0) 


0.141 


0.134 


-0.002 


0.006 


0.0 


0.3 


P ■ 


ud 1 35*1 


0.775(0.4) 


0.790 


0.781 


-0.015 


-0.005 


1.9 


0.2 


hi 


: ud 1 


1.230(3.2) 


1.283 


1.243 


-0.053 


-0.014 


2.5 


0.2 


Oi 


: 1 


1.230(40.) 


1.425 


1.320 


-0.195 


-0.090 


0.2 


0.1 


7T : 


itrf 2iS'o 


1.300(100) 


1.493 


1.435 


-0.193 


-0.135 


0.0 


0.0 


a2 


: 1 3P2 


1.318(0.6) 


1.276 


1.310 


0.042 


0.008 


12.9 


0.5 


P ■ 


ud 2 


1.465(25.) 


1.745 


1.684 


-0.280 


-0.219 


1.3 


0.8 


ao 


: ud 1 3Po 


1.474(19.) 


1.165 


1.024 


0.309 


0.450 


2.6 


5.6 




: ud 1 i_D2 


1.672(3.2) 


1.815 


1.763 


-0.143 


-0.090 


18.2 


7.2 


as 


: ud 1 ^Ds 


1.689(2.1) 


1.663 


1.718 


0.026 


-0.029 


1.2 


1.5 


ai 


: ud 1 ^Di 


1.720(20.) 


1.944 


1.847 


-0.224 


-0.127 


1.2 


0.4 


0.2 


: ud 2 3P2 


1.732(16.) 


2.025 


2.009 


-0.293 


-0.277 


3.3 


3.0 


TT : 


ud 3 ^.So 


1.816(14.) 


2.090 


2.037 


-0.274 


-0.221 


3.8 


2.5 


h2 


: ud 2 


1.895(16.) 


2.300 


2.267 


-0.405 


-0.372 


6.4 


5.4 


a4 


: ud 1 


2.011(12.) 


1.984 


2.057 


0.027 


-0.046 


0.1 


0.1 


62 


: ud 3ID2 


2.090(29.) 


2.704 


2.700 


-0.614 


-0.610 


4.5 


4.4 


P ■ 


ud 3 "'^i 


2.149(17.) 


2.281 


2.326 


-0.132 


-0.177 


0.6 


1.1 


03 


: ud 2 3£)3 


2.250(45.) 


2.275 


2.290 


-0.025 


-0.040 


0.0 


0.0 


as 


: ?id 1 '■'Gs 


2.330(35.0) 


2.258 


2.349 


0.072 


-0.019 


0.0 


0.0 


ae 


: (id 1 -'i/e 


2.450(130) 


2.500 


2.609 


-0.050 


-0.159 


0.0 


0.0 



Table 2: ud Mesons, Theory 1 and 2 - In this table and the ones below, the meson masses are in units of GeV, 
with experimental errors given parenthetically in units of MeV. 



SU, 


sd Mesons 


Exp. 


Thl. 


Th2. 


Exp.-Thl. 


Exp.-Th2. 


X'-Thl. 


X'-Th2. 


K- 


: sul'-So 


0.494(0.0) 


0.485 


0.519 


0.008 


-0.025 


0.7 


6.4 


K° 


: sdl^So 


0.498(0.0) 


0.488 


0.520 


0.010 


-0.022 


1.0 


5.0 


K* 


- :sul ^Si 


0.892(0.3) 


0.917 


0.896 


-0.025 


-0.004 


6.0 


0.2 


K* 


- : sd 1 ^Si 


0.896(0.3) 


0.919 


0.897 


-0.023 


-0.001 


4.8 


0.0 


K- 


: SU 1 iPi 


1.272(7.0) 


1.356 


1.339 


-0.084 


-0.067 


1.4 


0.9 


K* 


- : SU 1 3Pi 


1.403(7.0) 


1.419 


1.359 


-0.016 


0.044 


0.1 


0.4 


K* 


- : SU 23^1 


1.414(15.) 


1.759 


1.706 


-0.345 


-0.292 


5.3 


3.8 


K* 


" : SU 1 3Po 


1.425(50.) 


1.132 


1.079 


0.293 


0.346 


0.3 


0.5 


K* 


- : SU 1 3P2 


1.426(50.) 


1.379 


1.404 


0.047 


0.022 


6.8 


1.4 


K* 


- : sd 1 3P2 


1.432(1.3) 


1.380 


1.405 


0.052 


0.027 


10.2 


2.8 


K- 


: SU 2i5'o 


1.460(40.) 


1.523 


1.476 


-0.063 


-0.016 


0.0 


0.0 


K* 


- : SU 1 3p»i 


1.717(27.) 


1.922 


1.837 


-0.205 


-0.120 


0.6 


0.2 


K- 


: SU 1 ID2 


1.773(8.0) 


1.835 


1.803 


-0.062 


-0.030 


0.6 


0.1 


K* 


- : SU 1 '^Ds 


1.776(7.0) 


1.740 


1.792 


0.036 


-0.016 


0.3 


0.0 


K* 


- -.sul ^D2 


1.816(13.) 


1.824 


1.795 


-0.008 


0.021 


0.0 


0.0 


K- 


: SU 3 ^So 


1.830(13.) 


2.115 


2.081 


-0.285 


-0.251 


0.5 


0.4 


K* 


- : su2 3P2 


1.973(33.) 


2.078 


2.060 


-0.105 


-0.087 


0.1 


0.1 


K* 


- : SU 1 


2.045(9.0) 


2.045 


2.117 


0.000 


-0.072 


0.0 


0.6 


K* 


- : SU 2 ^D2 


2.247(17.) 


2.326 


2.313 


-0.079 


-0.066 


0.2 


0.1 


K* 


- : SU 2 


2.324(24.) 


2.642 


2.600 


-0.318 


-0.276 


1.7 


1.3 


K* 


- : SU 1 


2.382(14.) 


2.309 


2.401 


0.073 


-0.019 


0.3 


0.0 


K* 


- : SU 2 


2.490(20.) 


2.555 


2.600 


-0.065 


-0.110 


0.1 


0.3 



Table 3: su and sd Mesons, Theory 1 and 2 



10 



ss Mesons 


Exp. 


Thl. 


Th2. 


Exp.-Thl. 


Exp.-Th2. 


Y^-Thl. 


Y^-Th2. 

A. 


6- ss 1 


1 niQCn 01 


1.050 


1.013 


-0.031 


0.006 


9.3 


0.4 


(h • ss 1 3 Pn 


1 370fl00l 


1.211 


1.175 


0.159 


0.195 


0.0 


0.0 


(j) : ss l^Pi 


1.518(5.0) 


1.480 


1.437 


0.038 


0.081 


0.5 


2.5 


(f>: ss 1^P2 


1.525(5.0) 


1.496 


1.506 


0.029 


0.019 


0.3 


0.1 


6 ■ ss 2 


1 680f20 ) 


1.811 


1.875 


-0.131 


-0.195 


0.4 


0.9 


(?!) : ss 1 ^Ds 


1.854 7.0 


1.839 


1.879 


0.015 


-0.025 


0.0 


0.1 


(f>: ss2 3P2 


2.011(70) 


2.149 


2.128 


-0.138 


-0.117 


0.0 


0.0 


: ss 3 3P2 


2.297(28.) 


2.612 


2.603 


-0.315 


-0.306 


1.3 


1.2 



Table 4: ss Mesons, Theory 1 and 2 



CV.. ( 




Exp. 


Thl. 


Tli2. 


Exp.-Thl. 


Exp.-Th2. 


A --Thl. 


A^-Th2. 


: 


cu I^Sq 


1.865(0.2) 


1.865 


1.876 


0.000 


-0.011 


0.0 


1.1 


D+ : 


cd 1 ^So 


1.870(0.2) 


1.872 


1.883 


-0.003 


-0.013 


0.1 


1.7 




: cu 1 ^Si 


2.007(0.2) 


2.013 


2.007 


-0.006 


0.000 


0.4 


0.0 


D*+ 


: cd 1 ^Si 


2.010(0.2) 


2.019 


2.013 


-0.008 


-0.002 


0.7 


0.1 


D*° 


: cu 1 3Po 


2.352(50.) 


2.224 


2.221 


0.128 


0.131 


0.1 


0.1 


D*+ 


: cd 1 3Po 


2.403(14.) 


2.232 


2.230 


0.171 


0.173 


1.5 


1.5 


D+ : 


cd 1 3P2 


2.460(3.0) 


2.398 


2.414 


0.062 


0.046 


3.9 


2.1 


D*° 


: CM 1 3P2 


2.461(1.6) 


2.393 


2.409 


0.069 


0.052 


13.2 


7.7 



Table 5: cu and cd Mesons, Theory 1 and 2 



cs Mesons 


Exp. 


Thl. 


Th2. 


Exp.-Thl. 


Exp.-Th2. 


X^-Thl. 


X^-Th2. 


Z^s : cs 1 '^So 


1.968(0.3) 


1.972 


1.974 


-0.004 


-0.006 


0.1 


0.3 


D*:csl ^Si 


2.112(0.5) 


2.138 


2.119 


-0.026 


-0.007 


5.4 


0.4 


D*:csl 3Po 


2.318(0.6) 


2.348 


2.340 


-0.031 


-0.022 


6.9 


3.5 


pis : cs 1 iPi 


2.535(0.3) 


2.505 


2.499 


0.030 


0.036 


8.3 


11.6 


P; : cs 1 3P2 


2.573(0.9) 


2.534 


2.532 


0.039 


0.040 


8.4 


8.9 


P* : cs 2 


2.690(7.0) 


2.714 


2.702 


-0.024 


-0.012 


0.1 


0.0 




Table 6: cs 


Mesons, Theory 1 and 2 






cc Mesons 


Exp. 


Thl. 


Th2. 


Exp.-Thl. 


Exp.-Th2. 


X'-Thl. 


X'-Th2. 


77c : cc 1 


2.980(1.2) 


2.965 


2.973 


0.015 


0.007 


1.0 


0.2 


JliilS) : cc l^Si 


3.097(0.0) 


3.131 


3.128 


-0.034 


-0.031 


11.4 


9.7 


Xo : cc 1 3Po 


3.415(0.3) 


3.395 


3.397 


0.020 


0.018 


3.7 


3.0 


Xi : cc 1 3Pi 


3.511(0.1) 


3.506 


3.505 


0.005 


0.006 


0.2 


0.4 


hi-.cZl iPi 


3.526(0.3) 


3.522 


3.523 


0.004 


0.003 


0.2 


0.1 


X2 : cc 1 3P2 


3.556(0.1)) 


3.562 


3.557 


-0.006 


-0.001 


0.4 


0.0 


ri^:cc2 ^So 


3.637(4.0) 


3.606 


3.602 


0.031 


0.035 


0.6 


0.7 


tJj(2S) : cc2^Si 


3.686(0.0) 


3.688 


3.689 


-0.002 


-0.002 


0.0 


0.1 


^j{lD) : cc l^Di 


3.773(0.4) 


3.807 


3.807 


-0.034 


-0.034 


0.9 


0.9 


X2 :cc2 3p2 


3.929(5.0) 


3.980 


3.983 


-0.051 


-0.054 


1.0 


1.1 


V'(35) : cc3^Si 


4.039(10.) 


4.086 


4.092 


-0.047 


-0.053 


0.2 


0.3 


V'(2P) : cc 2^Di 


4.153(3.0) 


4.164 


4.169 


-0.011 


-0.016 


0.1 


0.3 


xp{AS) : cc 4^51 


4.421(4.0) 


4.410 


4.426 


0.011 


-0.005 


0.1 


0.0 


V'^SP) : cc S^Pi 


4.421(4.0) 


4.467 


4.483 


-0.046 


-0.062 


1.2 


2.3 


xJj{5S) : cc 5^51 


4.800(100) 


4.690 


4.719 


0.110 


0.081 


0.0 


0.0 


V'(4P) : cc 43Pi 


4.880(100) 


4.735 


4.764 


0.145 


0.116 


0.0 


0.0 


V'(65) : cc 6 35-1 


5.180(100) 


4.940 


4.983 


0.203 


0.197 


0.0 


0.0 


t/'(5P) : cc b^Di 


5.290(100) 


4.977 


5.020 


0.350 


0.270 


0.1 


0.1 



Table 7: cc Mesons, Theory 1 and 2 
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bu, 


bd bs Mesons 


Exp. 


Thl. 


Th2. 


Exp.-Thl. 


Exp.-Th2. 


X^-Thl. 


X^-Th2. 


B- 


: bul'So 


5.279(0.3) 


5.281 


5.283 


-0.002 


-0.004 


0.0 


0.2 


B° : 


bd 1 ^So 


5.280(0.3) 


5.282 


5.284 


-0.003 


-0.005 


0.1 


0.2 


B*- 


: bu 1 ^Si 


5.325(0.5) 


5.335 


5.333 


-0.010 


-0.008 


0.8 


0.5 


B*- 


: bu 1 


5.747(2.9) 


5.671 


5.687 


0.076 


0.059 


6.2 


3.8 




bs 1 ^Sq 


5.366(0.6) 


5.373 


5.367 


-0.007 


-0.001 


0.3 


0.0 


Bf 


: bs 1 ^Si 


5.413(1.3) 


5.441 


5.430 


-0.029 


-0.017 


3.0 


1.0 


Bf 


: bs 1 


5.829(0.7) 


5.789 


5.792 


0.040 


0.037 


10.9 


9.4 


Bf 


: 6s 1 


5.840(0.6) 


5.805 


5.805 


0.035 


0.035 


8.9 


9.0 


Be 


: be 1 


6.276(21.) 


6.249 


6.251 


0.027 


0.025 


0.4 


0.4 



Table 8: bu , bd and 6s Mesons, Theory 1 and 2 



66 Mesons 


Exp. 


Thl. 


Th2. 


Exp.-Thl. 


Exp.-Th2. 


X^-Thl. 


X^-Th2. 


■qi^-.bbl ^Sq 


9.389(4.0) 


9.337 


9.330 


0.052 


0.059 


1.6 


2.0 


T(15) : 66 l^^i 


9.460(0.3) 


9.444 


9.444 


0.016 


0.017 


2.4 


2.6 


Xw ■■ bb 1 


9.859(0.4) 


9.836 


9.834 


0.023 


0.026 


4.6 


5.6 


Xbi-bbl^Pi 


9.893(0.3) 


9.886 


9.886 


0.007 


0.007 


0.4 


0.4 


Xb2-bbl^P2 


9.912(0.3) 


9.922 


9.920 


-0.010 


-0.008 


0.9 


0.6 


T{2S) : bb 2^ Si 


10.023(0.3) 


10.022 


10.022 


0.001 


0.002 


0.0 


0.0 


T{D) : 66 2 3£>2 


10.161(0.6) 


10.178 


10.179 


-0.017 


-0.018 


2.1 


2.3 


Xbo-bb2^Po 


10.232(0.4) 


10.230 


10.229 


0.002 


0.003 


0.1 


0.1 


Xbi :66 2 3Pi 


10.255(0.5) 


10.261 


10.262 


-0.006 


-0.007 


0.3 


0.4 


Xb2-bb2^P2 


10.269(0.4) 


10.284 


10.286 


-0.015 


-0.017 


1.9 


2.5 


T{3S) : 66 3 ^S"! 


10.355(0.6) 


10.366 


10.368 


-0.011 


-0.013 


0.8 


1.2 


T{4S) : bbi^Si 


10.579(1.2) 


10.626 


10.633 


-0.046 


-0.053 


8.8 


11.7 


T{5S) : bb5^Si 


10.865(8.0) 


10.844 


10.857 


0.021 


0.008 


0.1 


0.0 


T{6S) : 66 6 351 


11.019(8.0) 


11.036 


11.055 


-0.017 


-0.036 


0.0 


0.2 



Table 9: 66 Mesons, Theory 1 and 2 



Meson Family 






# Mesons 


Average x^ 


Thl 


Th2 


Thl 


Th2 


ud 


62.1 


30.5 


20 


3.1 


1.5 


su , sd 


37.7 


26.5 


22 


1.8 


1.3 


ss 


11.0 


4.5 


8 


1.4 


0.6 


cu , cd 


20.9 


14.6 


8 


2.6 


1.8 


cs 


29.4 


26.6 


6 


4.9 


4.4 


cc 


28.2 


23.4 


18 


1.6 


1.3 


bu , bs , be 


31.2 


25.9 


9 


3.5 


2.9 


bb 


24.8 


25.1 


14 


1.8 


1.8 


Total 


245.3 


177.1 


105 


2.4 


1.7 



Table 10: x^ by Family for Thl and Th2 
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1. CLiiiiiV 


Exp. 


Thl. 


T1l2. 


ud 


DOO 


04y 


04 < 


su 


OijO 




■^77 


oci 


oyo 




■^77 


cu 




1 4S 


lOl 


Cd 


1 /in 




1 "in 

loU 


cs 


1/1/1 

144 


100 


1 /I ^ 

140 


cc 


117 


100 


1 

100 


bu 


40 


f;4 


ou 


s's 


47 


68 


63 


bb 


71 


107 


114 


Table 11: 


Ground State Sinj 


^let- Triplet Splitting 


5S (MeV) 



shifts from about a factor of 2 to almost even. The heaviest mesons are also the smallest in mean radius. This 
means that they are less likely to experience the effects of the S'^ and —V'^ portions of the confining interactions. 
As the mesons become large, they experience more of the effects of these parts of the potentials. The most 
dramatic improvement from the inclusion of the timelike vector confining potential V{r) (Th2) is with the light 
quark ud family. Referring now to Tables most of the improvement comes from that of the fits to the 02, 62 
and bi mesons. For the sm, sd family the largest improvement comes from the lowest lying '^P2 mesons. The 
ground state singlet-triplet splitting changes from an overestimation to an underestimation. For the ss family 
the most significant improvement is in the ground state, although somewhat off-balanced by a worse fit for the 
lowest lying '^Pi state. In the case of the cu, cd family the main improvement is from the '^P2 mesons. For 
the as mesons there is a slight overall improvement for Th2 with offsetting changes for the ^Pq and ^Pi mesons. 
There is only a very slight improvement for the cc mesons. For both Thl and Th2 the worst fit is to the J/ip 
meson with a mass too large by about 30 MeV. It cannot be adjusted downward by lowering the charm mass 
due to the fact that other mesons in this family would be pushed further from the data. With the heavy-light 
family, a single b quark, there again is not much overall change and even less in the 66 family although another 
significant improvement is in the '^P2 bu state. An oddity with the 66 is the sudden increase in at the 4^5i 
meson, the worst fit of all the mesons in terms of the incremental x^- Since that meson is closest to threshold, 
its mass will be most affected by it, whereas our theoretical model does not take threshold effects into account. 

A possible explanation of why most improvements come for the states is that the effect on the spin-orbit 
coupling due to the Thomas terms is opposite in sign for the timelike vector and scalar mesons. Without the 
balancing effect of the timelike vector confining interaction, the scalar interaction enforces an inversion of the 
spin-orbit splittings of the light mesons that are far too distorted for the ud and su multipletj^. Also, the 
long range scalar parts contribute oppositely in sign from the short range vector part attributable to the A[r) 
potential. 

We now examine another important feature of our method: the goodness with which our equations account 
for spin-dependent effects (both fine- and hyperfine- splittings). Table [TT] shows the best fit versus experimental 
ground state singlet-triplet splittings and six versus four of the ten hyperfine splittings are improved using Th2 
over Thl. Both give good fits for all hyperfine ground state splittings except for the rj^ — ip system and rj^, — T 
system which for the latter over estimate the splittings by about 50%{3 One problem with the fit for the cc 
system of mesons may be due to the fact that the D* ^P2, ^Pi and D* ^P2 fits are significantly low while the 
J /ip fit is significantly high. Lowering the c quark mass corrects the J /ip mass while raising the £>*, D* P state 
masses would require raising the c quark mass. Reducing one discrepancy would worsen the other, at least in 
our three invariant function approach 

For the spin-orbit splittings Table [T^ gives the R ratios ('^P2 — ^ Pi)/(^Pi — ^ Po))- Both sets of fits are very 
poor for the two lightest multiplets. The fact that Thl has the same sign for the ud as the experiment values is 
not an indication that it gives reasonable results since the negative sign originates from the numerator instead of 
the denominator. Of the four remaining multiplets, Th2 gives a better fit on 3. It must be said, however, that 

Another possible source of the strange multiplet inversion, is that the observed ^Pq states of 1450 for the ud and 1430 for the 
us systems are, in fact not zero node states, but rather one node excited states. This may, in our formalism, give room to an 
interpretation of the uJ(980) and the k(700 — 900) mesons as possible candidates for the zero node states. Our tables support that 
more than the identification of a zero node 1474 for the ud and 1425 for the us systems. However, using the parameters on our 
model, we would obtain the 2^Po values of 1800 and 1850 for those one node states, well above the 1474 and 1425 experimental 
values. So this docs not appear to be a plausible alternative for either Thl or Th2. 

^''We point out, however, that our model does much better in simultaneously working with heavy and light qq hyperfine splittings 
than that obtained by typical constituent (non-chiral) quark models (a major exception discussed here in detail is the quasipotential 
model of [Tj) 



13 



1. Ojiiiii V 


Exp. 


Thl. 


Th2. 


ud 


-U.oD 


-U.O ( 


-U.Uo 


SU 


1 

-i.UO 


n 1 /I 


U.iD 


SS 


u.uo 


u.uu 


0.26 


cc 


n 47 


n "lO 


0.' 


18 


bb 


(l'^2,i.o) 


0.56 


0.72 


0.65 


bb 


(2'^2,i,o) 


0.61 


0.74 


0.73 


Table 12: 


Spin-Orbit Splitting R Ratios 








Family 


Exp. 


Thl. 


Tli2. 






ud 


76 


30 


36 






su 


146 


9 


14 






cc 


-1 


3 


-1 





Table 13: Splitting Between ^Pi and weighted triplet states (MeV) 



none of the better fits are very good except for the cc. From the experimental point of view the poor R value for 
the ud and us may be the uncertain status of the ^Pq light quark meson bound states, or theoretically our low 
^Po theoretical meson masses. Also, the lack of any mechanism in our model to account for the effects of decay 
rates on level shifts undoubtedly has an effect. Another likely cause for the poorer performance of Thl as one 
goes from heavier to light mesons is that the radial size of the meson grows so that the long distance interactions, 
in which the scalar interaction becomes dominant, play a more important role. This effect is blunted by Th2 
as seen in the ud numbers. In Table [2] The spin-orbit terms due to scalar interactions are opposite in sign and 
tend (at long distance) to dominate the spin-orbit terms due to vector interactions for Thl, but less so in Th2. 
This results in partial to full multiplet inversions as we proceed from the ss to the ud mesons. This inversion 
mechanism is less for Th2 than for Thl because of the value of 

The hyperfine structure of our equations also influences the splitting between the ^Pi and the weighted sum 
[5(^P2) + 3('^Pi) + l('^Po)]/9 of bound states. Table [T3l indicates the agreement of the theoretical and experimental 
mass differences is excellent for the cc system, too small but of the right sign for the us system and ud systems. 
The agreement, however, for the light systems is nevertheless considerably better than that in the case of the 
fine structure splitting R ratios. Note that in the case of unequal mass P states, our calculations of the two 
values incorporate the effects of L • (cri — (T2) and L • cri x cr2 which mix spin. 

Next consider the mixing due to the tensor term between orbital D and radial S excitations of the spin-triplet 
ground states. This mixing occurs most notably in the cc,us and ud systems. Table [HI show that Thl is better 
than Th2 although both are pretty far off the mark. For the charmonium system, the lower doublet results are 
high whereas the higher doublet results are low. 

Next we consider the effects of the change from Thl to Th 2 on the radial excitations. Table [12] shows that 
for the most part Thl gives slightly better results for the radially excited states, although where both theories 
are furthest off (the ud states) Th2 gives better results. The radially excited ud mesons have a larger mean 
radius than for the heavier meson and thus the temporizing effects of the —V^ term tends to counteract more 
the increased confining potential for large r from linear to quadratic due to the terms. 

Finally we comment on the isospin splittings shown in table 1161 There are two effects we must consider here: 
the positive d — u mass differences of about 6 MeV for both theories and the Coulomb interaction between the 
quarks on the order of a x 197 MeV or less depending on the meson sizes. The Coulomb interaction is counter 
to the d — u mass difference for the sd — su and bd — bu splittings while enhancing the d — u mass difference for 
the cd — cu splittings. These alternatively competing and enhancing effects are seen in the sizes of the splittings 



Family 


Exp. 


Thl. 


Th2. 


ud 


255 


199 


163 


su 


303 


163 


131 


cc il^Di-2^Si) 


87 


119 


118 


cc {2^Di-3^Si) 


114 


78 


77 



Table 14: Mixing Due to the Tensor Term Between Orbital D and Radial S Excitations of the Spin- Triplet 
Ground States (in MeV) 
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SI n 
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^ ■ 
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/ Dl 




cc 


Vc ■ 


1 9I i?„ 
i, Z Oq 


00 / 




fi9Q 
uzy 


cc 


ih ■ 

'r • 


1 2 3^51 


589, 353 


557, 398 


561, 403 


bb 


T : 


1,2, 3, 4, 5,6351 


563, 332, 224, 286, 154 


578, 344, 260, 218, 192 


578, 346, 265, 224,198 



Tabic 15: Radial Excitations (McV) 



Family 




Exp 


Theory 1 


Theory 2 


sd — su : 




4 


3 





sd — su : 




4 


2 


1 


sd — su : 




6 


1 


1 


cd ^ cu : 


I'So 


5 


7 


7 


cd — cu : 


13^1 


3 


6 


6 


cd — cu : 


l3Po 


51 


8 


9 


cd — cu : 


1=^P2 


1 


-5 


-5 


bd — bu : 


I'So 


1 


2 


1 


d — u mass 




5.8 


6.0 



Tabic 16: Isospin Splitting (McV) 



for both theories as you read down the table from the sd — su through the cd — cu to the bd — bu splitting. 
For the K — K* family the values for the isospin splittings are 3 and MeV for Thl and 2 vs the experimental 
value 4 MeV for the singlet ground states while for the triplet the isospin splittings are 2 and 1 MeV vs the 
experimental value 4 MeV. The experimental splitting grows for the orbital excitation {K2) to 6 MeV. The 
probable reason for the increase is that at the larger distances, the influence of the Coulomb differences becomes 
small so that only the d — u mass difference influences the result. Our theories do not show a similar increase for 
the orbital excitations. In the case of the — 1)° splitting our mass differences for Thl and Th2 are 7 and 7 
MeV respectively versus the experimental mass difference of just 5 MeV. Here we see the opposite overall effect 
between the combined effects of the d — u mass difference and the slightly increased electromagnetic binding 
present in the case of the 1?° and the slightly decreased binding in the case of the D'^. Whereas in the kaon 
system the results are too small, for the D the results are too large . This can be partially understood since 
the Coulomb and d — u mass differences work in concert with the Coulomb potential for these doublets. These 
effects work in the same way for the spin-triplet splitting resulting in the theoretical values of 6 and 6 MeV 
for the two theories compared with the experimental value 3 MeV. For the isodoublet we obtain -5 and -5 
MeV versus about 1 for the experimental value again showing the expected opposite trend from that of the kaon 
system. The experimental splitting between the '^Pq isodoublet of 51 MeV appears incomprehensibly large. Our 
two values are 8 and 9 MeV. The isospin splittings that we obtain for the spin singlet B meson system are 2 
and 1 MeV for Thl and Th2 versus 1 MeV. Here the competing effects cancel as in the kaon system only more 
so since the mesons are smaller and thus the Coulomb parts play a stronger role than for the kaon. 

5 The Effective Relativistic Schrodinger Equation with Flavor Mix- 
ing for Spin-Zero Isoscalar Mesons. 



Consider the general eigenvalue equation p5p for an isoscalar meson, one with quark structure qq. As seen in 
Appendix A the mass dependence appearing in directly or indirectly through m^, £„,, ei, £2, is of four types: 



TOl77l2, TOj 

quark- antiquark pairs, 
shorthand as 



and 



The actual isoscalar mesons consist of mixtures of three equal mass 



We write the three separate equal mass versions of ((35l) . using Eq. pO|) . together in 



[p^ + $^(r,TOi = 7712, w,cri, 0-2)] 



-{u? -AM?) 



i'dd 



(37) 
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Parameter Thl Th2 

~lrr^, OlOCM 0.1070 

6md 0.1378 0.1055 

Sm, 0.0468 0.0578 

m,(^o) 0.0737 0.1015 

m,(^) 0.2175 0.2261 

m^fa,) 0.4297 0.4536 

Table 17: Mixing Parameters for Th 1 and 2 (GeV) 



Mesons Exp. Thl. Th2. Exp.-Thl. Exp.-Th2. x^-Thl. x^-Th2. 

tt'J : 1 0.135(0.0) (U39 0134 41004 (0)01 (L2 OlT" 

r;: l^S-Q 0.548(0.0) 0.548 0.548 0.000 0.000 0.0 0.0 

r]' : 1 ^Sq 0.958(0.2) 0.958 0.958 0.000 0.000 0.0 0.0 



Table 18: cs Mesons,Theory 1 and 2 (GeV) 



in which 


















rrid 
















(38) 



Eq. p7p does not include mixing between the pairs. Motivated by ideas presented by Brayshaw[^ [5, we model 
the effects of qiQi — > qjcjj via two gluon annihilation and creation as an effective scalar potential by postulating 
a symmetric matrix M that is not diagonal, 



ruu 





" 




■niu 





" 




Srriu 





rud 





^ M = 





md 





+ 


^SmuSnid 








ms_ 




_ 









y/SniuSms 



y/SmdSnis 



\/SmdSms 
Srus 



(39) 



Suppose that an orthogonal matrix M diagonalizes 



Then Eq. ^ becomes 



16 



(40) 



$^(r,Mz))]7/; = -(u;2 



4M^)V'. 



This gives us, in essence, three new effective families of equal quark-anti-quark mesons, like ones that contain 
&, c, s, M, d except that mixtures are involved. In this paper we see if this idea is successful for the ground state 
pseudoscalar isoscalar family of mesons alone. With the three parameters, one obtains three different effective 
quark masses, one for each isoscalar family. The three Srrii are adjusted to give the best fit to the correct 
7r°, 77, rj' masses. ^ Table [ni gives the values of 5mu, Snid, 6ms together with the three effective quark masses, 
the eigenvalues of M which wc call m^j-^o), 'm.gj^jyy 



Paralleling the earlier Tables [2]|9l Table [TSl gives the best fit values for the 7r°, 77, and 77' mesons. The 
predicted quark content becomes a further test of our model. 



^^Brayshaw considered the modification of the meson mass by wl - 
subspace. This is equivalent to the mixing matrix 



wl + \G)Z{G\ with Z a fixed parameter and \G) in the ns 



T =Z 

with the property that (for real matrix elements) 



{nn\G){G\nn) {nn\G){G\ss) 
{ss\G){G\nn) {sS\G){G\sS) 





til 


tl2 




*21 


t22 



illt22 = 412*21 = 1*12!^ 



^12> 



just as with the choice of Ea.l l39l l. 

'^^We ignore here the coupling that would result from this diagonalization that would be brought on by the Coulomb interactions 
between the equal mass qq pairs. 

^'^ These fits are simultaneous with the fits of the earlier 105 mesons, with the same quark masses and potential parameters used in 
the generation of Tables II2I9I I. Oddly, a precise fit to the tt" was not possible in either theory, in spite of the three extra parameters 
Srriu, Sm^, and Srris available. 
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The corresponding eigenvectors are quite close to the mixtures 



W) = 



1 

cosf 
sin 



1 
1 

-2 

1 
1 

-2 



17^3), 



surf 



V3 



cosy 



cosy |?7g) — smy \rji) 



sin0|?7g) + COS0 |?7i), 



(41) 



With the three eigenvectors in matrix form we find 





' 0.770 


0.470 





43f 


\v) W)] = 


-0.638 


0.574 





514 




-0.006 


-0.671 





742 



(42) 



corresponding to 6 = —12.6 degrees for Thl and 





' 0.717 


0.542 





438" 


\v) W)] = 


-0.697 


0.565 





442 




-0.008 


-0.622 





783 



corresponding io 9 — —16.3 degrees for Th2. 
corresponding to using the formula [44] . 

tan^e 



(43) 



The Th2 value is consistent with chiral perturbation theory results 



[quad] 



— 3to^ 



3m^, 



(44) 



With the meson masses listed in Tables [5] and [3] for Th2 we obtain 9\^q^g^j^ — —17.2 degrees, reasonably close to our 
Th2 result of —16.3 degrees. On the other hand, using the values listed there for Thl gives (?[g„cic;] = —7.3 degrees 
which is significantly different from —12.6 degrees. Using the experimental masses gives ^[q„od] — —11.5 degrees. 
Radiative vector meson decays give an angle between —10 and —20 degree while fits to tensor decay widths give 
— 17 degrees. Note that even though the chiral symmetry and its breaking are not built into our model as in 
for example it is in [35] , it is of interest that the use of Eq. pij) does produce a consistent result by use of our 
theoretically computed masseJ^. 



6 The Quasipotential Equation of Ebert, Faustov and Galkin 
6.1 The Model and Comparisons with TBDE 

The model which we critically examine here gives excellent fits to the meson spectrum as well as numerous meson 
decay rates. The quasipotential approach of Ebert, Faustov and Galkin (EFG) [7] is a local one very similar to 
that of Todorov [TU] and Aneva, Karchev, and Rizov [3B], discussed in Sec. 2 (see Eq. ([^ l. Insofar as our 
discussion given in that section is concerned, the main difference between the TBDE and EFG approach is the 
replacement of the timelike vector confining interaction in Eq. (j26p with a different confining vector interaction: 



v(p-k)/3i/32 ^ v(p-k)ri,,r^. 

They include as do we a scalar confining potential. In coordinate space their choice is 

V{r) = {l-e){Ar + B), 

S{r) = e{Ar + B). (46) 



^^There are at least three other additional hints of a possible emergent chiral symmetry and its breaking in the TBDE model . 
a) the relatively small quark u, d, masses ~ 60 — 100 MeV compared with ~ 300 MeV with most other quark models b) Our tt mass 
Bases tow 
i[39],[l3] 
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For their electromagnetic-like vector interaction they use the Coulomb gauge (instead of the Feynman gauge 
used in the TBDE). 

4 

Its momentum space form is 

(P - k) {p~ k) (p - k) 

The addition of the Pauli term with their value k = — 1 has the effect of cancelling the lowest order spin-dependent 
contributions in each factor of F^^ when sandwiched between on energy shell spinors. In the nonrelativistic limit 
their scalar and vector confining interactions combine to {Ar-\-B) with A = 0.18 GeV^, B=— O.SGeV. In essence 
their approach embodies a modified version of the Cornell potential into the local quasipotential approach. Their 
choice of k was fixed by an analysis of the fine structure splitting of heavy quarkonia ^Pj states I 46j and their 
choice of e = — 1 is determined from considerations of charmonium radiative decay. 



We compare in Tables [12] and HOI their results to our Th 2 by listing the deviations from the experimental 
results and respective . Of the 86 common mesons fit to the respective models, the collected results of 
quasipotential approach of [7], [1B]-[5D] (EFG) are more accurate in 69 of the fits (61 when comparing x^0- 
This includes the difficult singlet-triplet splittings for the ground and excited states of charmonium and the 
ground state of bottomonium, as well as good fits to most of the radial and orbital excitations of the ground 
states of the light mesons. Particularly interesting examples for the tt — p family are the ^Pi Pi splitting and 
the radial excitation of the singlet and triplet S*— states, and for the K-K* families the ^Pi — ^ Pi splitting and 
the '^Po mass. These three areas of their spectrum are noteworthy improvements over the TBDE approach. It 
is hard, however, to give an even theoretical comparison between the two approaches for a number of different 
reasons. It is worthwhile, however, to point out the differences in the two approaches, summarized in Table [21] 
First of all in our approach we give an overall fit to the entire spectrum. The approach of EFG to the spectrum 
is spread over several papers and it is not clear that a uniform parametrization would yield the same results as 
given in the tables (summarized here) from their separate papers. They do use the same values for the constants 
A, B, K, and e as well as the quark masses in the various papers. However, the A parameter they use in the 
parametrization of the short distance QCD-Coulomb part of the potential is different Another difference 
is that the static QCD potential in the Adler-Piran model displays explicitly the asymptotic freedom behavior 
by its radial dependence whereas in the Coulomb potential used by EFG, the asymptotic freedom behavior 
is displayed indirectly in their as by the quark mass dependence as seen in footnote 20. On the other hand, 
for their heavy quark bound states a radial modification displaying short distance QCD asymptotic freedom 
corrections was used but was not for the light quark bound states. Both the quasipotential approach of EFG 
and the TBDE used here, extending earlier work of Crater and Van Alstine, have three invariant interaction 
functions. Both use electromagnetic-like four-vector interactions; in the Feynman gauge for the TBDE and the 
Coulomb gauge for the quasipotential approach. In the EFG quasipotential approach the third interaction is a 
Pauli-modified vector interaction whereas in our approach it is timelike vector interaction. The potentials used 
in each of the three parts are of course different in the two approaches. The spin dependence, although similar 
for the most part have distinctly different origins. In the 16 component TBDE the kinematics and dynamics 
and spinors are tied together in one wave equation (see. e.g. Eq. P^ ). Its off shell dependence is fixed by the 
wave equation and the spinors are all interacting. In the quasipotential approach the potential is constructed 
in part from the actions of free particle spinors. This leaves substantial leeway in how the off-shell behavior is 
fixed. 

Both approaches have exact relativistic kinematics (from the use of 6^(w)) and do not use cither v/c or 
l/mq expansions. Both approaches lead to non-linear eigenvalue equations and give good values of the pion and 
kaon masses. Both approaches avoid singular effective potentials that would otherwise prevent nonperturbative 
spectral calculations. In the quasipotential approach of EFG, those singularities are avoided in an adhoc though 

^^Due to their large descrepancy on the (controversial) ^Pq -Di,(2370) meson, their of 389 is actually larger than our 255 (our 
fit here is altered to include less mesons and m„ = m^j so there are slight differences from the previous tables). If one eliminates 
that meson from the fit our reduces to 251 while their reduces to 69. It should be pointed out that their fit did not appear 
to be a least x^ one like our fit. It also was not an overall fits like ours. 

■^"^In their ealier papers where the fits to the heavier mesons are given they use = iir /{Pq ln(/i^/A^)) where ^ is the reduced 
mass and in the recent papers where the fits to lighter mesons are given they use as = 47r/(/3Q ln((/i'^ -|- Mg)/A-^)). In the former 
papers they use A = 169 or 178 MeV and in the recent ones A = 413 MeV. It is not clear how using just one form for all the mesons 
would affect the overall fit. 

^^In 1331 the Adler-Piran potential was replaced by a form that displays asymptotic freedom in the QCD coupling via 
(87r/27) /ln(K + B/{Ar)^). Although the model used there does not give as good a fit to the meson spectrum as the Adler- 
Piran model it does display asymptototic freedom in a simpler form. 
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Old Mesons 


Exp. 


TBDE 


EFG 


Exp.- TBDE 


Exp. - EFG. 


X'(TBDE) 


X^(EFG). 


TT : ud 1 ^Sq 


0.140 


0.134 


0.154 


0.005 


-0.014 


0.4 


2.5 


p : ud l^Si 


0.775 


0.779 


0.776 


-0.003 


-0.001 


0.1 


0.0 


bi-.udl 


1.230 


1.237 


1.258 


-0.007 


-0.028 


0.1 


0.9 


ai :ud 1 ^Pi 


1.230 


1.311 


1.254 


-0.081 


-0.024 


0.0 


0.0 


mud 2 ^So 


1.300 


1.426 


1.292 


-0.126 


0.008 


0.0 


0.0 


a2 -.ud 1 ^P2 


1.318 


1.303 


1.317 


0.015 


0.001 


2.0 


0.0 


p : ud2^Si 


1.465 


1.674 


1.486 


-0.209 


-0.021 


0.8 


0.0 


ao : ud I ^Pq 


1.474 


1.015 


1.176 


0.459 


0.298 


7.0 


3.0 


&2 : wrf 1 


1.672 


1.752 


1.643 


-0.080 


0.029 


6.8 


0.9 


as -.ud 1 ^-Ds 


1.689 


1.706 


1.714 


-0.017 


-0.025 


0.7 


1.4 


ai :ud 1 ^Di 


1.720 


1.836 


1.742 


-0.116 


-0.022 


0.4 


0.0 


02 : ltd 2 3P2 


1.732 


1.997 


1.779 


-0.265 


-0.047 


3.3 


0.1 


TT : ud 3 ^Sq 


1.816 


2.022 


1.788 


-0.206 


0.028 


2.6 


0.0 


&2 : 2 ID2 


1.895 


2.252 


1.960 


-0.357 


-0.065 


6.0 


0.2 


a4 : ud 1 ^F^ 


2.011 


2.042 


2.018 


-0.031 


-0.007 


0.1 


0.0 


&2 : wrf 3 iL»2 


2.090 


2.682 


2.216 


-0.592 


-0.126 


5.0 


0.2 


p : a^s"! 


2.149 


2.309 


1.921 


-0.160 


0.228 


1.1 


2.2 


uq : ud I ^Hq 


2.450 


2.590 


2.475 


-0.140 


-0.025 


0.0 


0.0 


K~ : sul ^So 


0.494 


0.528 


0.482 


-0.034 


0.012 


14.2 


1.6 


K*- : sul ^Si 


0.892 


0.898 


0.897 


-0.007 


-0.005 


0.5 


0.3 


K- -.sul^Pi 


1.272 


1.336 


1.294 


-0.064 


-0.022 


1.0 


0.1 


K*- : sul 


1.403 


1.354 


1.412 


0.049 


-0.009 


0.6 


0.0 


K*- : su2 ^Si 


1.414 


1.698 


1.675 


-0.284 


-0.261 


4.3 


3.6 


K*- -.sul 


1.425 


1.075 


1.362 


0.350 


0.063 


0.6 


0.0 


if* - : su 1 3P2 


1.426 


1.401 


1.424 


0.025 


0.002 


0.0 


0.0 


K- -.3112^80 


1.460 


1.414 


1.538 


0.046 


-0.078 


0.0 


0.0 


K*- :sul ^Di 


1.717 


1.828 


1.699 


-0.111 


0.018 


0.2 


0.0 


K- ■.sul^D2 


1.773 


1.795 


1.709 


-0.022 


0.064 


0.1 


0.8 


K*- -.sul ^Ds 


1.776 


1.784 


1.789 


-0.008 


-0.013 


0.0 


0.0 


K*- :sul 3£)2 


1.816 


1.787 


1.824 


0.029 


-0.008 


0.1 


0.0 


K ~ : su 3^So 


1.830 


2.069 


2.065 


-0.239 


-0.235 


4.1 


3.9 


K*- : su2^P2 


1.973 


2.050 


1.896 


-0.077 


0.077 


0.1 


0.1 


K*- : sul ^Fi 


2.045 


2.106 


2.096 


-0.061 


-0.051 


0.5 


0.4 


K*- : su2^D2 


2.247 


2.301 


2.163 


-0.054 


0.084 


0.1 


0.3 


K*- : su2 ^Fs 


2.324 


2.585 


2.348 


-0.261 


-0.024 


1.4 


0.0 


K*- -.sul ^Gs 


2.382 


2.387 


2.356 


-0.005 


0.026 


0.0 


0.0 


K*- : su2 ^Fi 


2.490 


2.585 


2.436 


-0.095 


0.054 


0.3 


0.1 


4>:ssl '^Si 


1.019 


1.017 


1.038 


0.002 


-0.019 


0.1 


4.1 


4>: ssl 


1.370 


1.175 


1.420 


0.195 


-0.050 


0.0 


0.0 


(j): ssl 3Pi 


1.518 


1.436 


1.464 


0.082 


0.054 


3.1 


1.4 


^ : ss 1 3P2 


1.525 


1.505 


1.529 


0.020 


-0.004 


0.2 


0.0 


(j): ss2 ^Si 


1.680 


1.868 


1.698 


-0.188 


-0.018 


1.1 


0.0 


(j): ssl 


1.854 


1.874 


1.950 


-0.020 


-0.096 


0.1 


2.2 


(j): ss2 3P2 


2.011 


2.120 


2.030 


-0.109 


-0.019 


0.0 


0.0 


(j): ssS 3P2 


2.297 


2.590 


2.412 


-0.293 


-0.115 


1.3 


0.2 



Tabic 19: Comparison of Th2{A,S,V) Fits with Ebcrt et al [7], [46]-[50] 



19 



New Mesons 


Exp. 


TBDE. 




Exp. - IBDE 


Exp.-EFG 


X -II3DE 


X -EIG 


D+ : cd 1 ^5*0 


1.870 


1.881 


1.871 


-0.012 


-0.001 


1.5 


0.0 


: cd 1 Di 


2.010 


2.010 


2.010 


A AAA 

0.000 


A AAA 

0.000 


A A 

0.0 


A A 

0.0 


D*+ : cd 1 ^Po 


2.403 


2.224 


2.406 


0.179 


-0.003 


2.0 


0.0 


D*° : cu 1 


2.460 


2.408 


2.460 


0.052 


0.000 


3.3 


0.0 


Ds-.csl ifi'o 


1.968 


1.976 


1.969 


-0.007 


-0.001 


0.5 


0.0 


D* : cs 1 


2.112 


2.120 


2.111 


-0.008 


0.001 


0.7 


0.0 


: cs 1 


2.318 


2.338 


2.509 


-0.020 


-0.191 


3.7 


320 


Ds-.c-sl 


2.535 


2.498 


2.536 


0.038 


-0.001 


15.4 


0.0 


D* : cs 1 


2.573 


2.531 


2.571 


0.042 


0.002 


11.8 


0.0 


P>* : cs 2 35-1 


2.690 


2.698 


2.731 


-0.008 


-0.041 


0.0 


0.4 


1 i~i 

T?c : cc 1 'So 


2.980 


2.972 


2.978 


0.008 


0.002 


0.3 


0.0 


J/^{IS) : cc 1 ^Si 


3.097 


3.126 


3.097 


-0.029 


0.000 


10.3 


0.0 


Xo-ccl 3Po 


3.415 


3.393 


3.423 


0.022 


-0.008 


5.2 


0.7 


Xi : CC 1 ■'Pi 


3.511 


3.500 


3.509 


0.010 


0.002 


1.3 


0.0 


/ii : cc 1 'Pi 


3.526 


3.519 


3.525 


0.007 


0.001 


0.6 


0.0 


X2 : CC 1 ^P2 


3.556 


3.553 


3.556 


0.004 


0.000 


0.1 


0.0 


Vc ■■ cc2'So 


3.637 


3.597 


3.663 


0.040 


-0.026 


1.1 


0.5 


tp{2S) : cc 2^ Si 


3.686 


3.683 


3.684 


0.004 


0.002 


0.1 


0.1 


V'(l-D) : CC l^Di 


3.773 


3.801 


3.795 


-0.028 


-0.022 


8.4 


5.2 


X2 :cc2 3p2 


3.929 


3.975 


3.972 


-0.046 


-0.043 


1.0 


0.9 


V'(3S') : cc S^S"! 


4.039 


4.083 


4.088 


-0.044 


-0.049 


0.2 


0.3 


yj{2U) : cc I L>i 


4.153 


4.160 


A -i Ci A 

4.194 


A AAT 

-0.007 


A A /I 1 

-0.041 


A 1 

0.1 


A 

2.0 


: bu 1 


5.279 


5.285 


5.280 


-0.006 


-0.001 


0.4 


0.0 


B^' : bd I '^So 


0.000 


0.000 


0.000 


0.000 


0.000 


0.0 


0.0 


D* — . U-TT 1 3 

is . bu I bi 


o.o25 


5.oo4 


O.OZD 


-0.009 


-0.001 


0.0 


A A 
0.0 




0. ( 4 ( 


iD.DO / 


^ 7/1 1 


U.UDU 


U.UUO 


4. 1 


U.U 


i)^ . OS i Do 


O.ODD 


0.0 <u 


0.0 ( Z 


n nn/i 

-U.UU4 


-U.UUu 


n 1 


U.o 


-Dg . ys 1 Ol 


A^'i 
0.^10 






m Q 


n 001 


1 7 


n 

U.U 


rj*0 • 1 3 p 


5.829 


5.793 


5.831 


0.037 


-0.002 


11.0 


0.0 


P*° : 6s 1 3P2 


5.840 


5.805 


5.842 


0.034 


-0.002 


10.5 


0.0 




Q '^89 




9 400 


055 


-0 01 1 

\J .yj ± ± 


2.1 


1 


TH .S""! • /)/} 1 3,^1 


9.460 


9.447 


9.460 


0.014 


0.000 


2.1 


0.0 




9.859 


9.835 


9.864 


0.024 


-0.005 


6.1 


0.2 


Y>,i : 66 1 3Pi 


9.893 


9.887 


9.892 


0.006 


0.001 


0.4 


0.0 




9.912 


9.921 


9.912 


-0.009 


0.000 


0.8 


0.0 


T(25) : 66 2 351 


10.023 


10.021 


10.020 


0.002 


0.003 


0.0 


0.1 


T(1D) : bbl^D2 


10.161 


10.178 


10.157 


-0.017 


0.004 


2.4 


0.1 


Xbo-bb2^Po 


10.232 


10.228 


10.232 


0.005 


0.001 


0.2 


0.0 


Xbi :66 2 3Pi 


10.255 


10.261 


10.253 


-0.005 


0.002 


0.3 


0.1 


XJ,2:&^'2 3P2 


10.269 


10.284 


10.267 


-0.015 


0.002 


2.4 


0.0 


T(3S') : 66 3 35*1 


10.355 


10.366 


10.355 


-0.011 


0.000 


1.0 


0.0 


T(4S') : 66435*1 


10.579 


10.628 


10.604 


-0.049 


-0.025 


11.7 


3.0 



Table 20: Comparison of Th2{A,S,V) Fits with Ebert et al [7], [46]-[50] 
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Properties 


TBDE 


Quasipotcntial Approach of EEC 


Invariant Interactions 


3-A(EM-like vector), S'(scalar), 
y (timelike vector) 


3- A(EM- like vector), S'(scalar), 
l^(Pauli-modified vector) 


QCD Coupling 


Coordinate Space Dependent 


Quark Mass Dependent 


Meson Fits 


Overall Spectrum, Same 
Parametrizations 


Spectrum in Parts with Different 
Parametrizations 


Spin- Dependence 


Fixed by the TBDE given A, S, V 


Fixed by the Quasipotcntial and A, S, V 


Kinematics 


Exact 


Exact 


Singular potentials 


Avoided by Dirac Equation Formalism 


Avoided by Adhoc Substitutions 


Numerical evaluations 


Yes 


Not for Heavy Mesons 


Chiral Symmetry 


Zero Pion Mass for niq — > 


Not Tested 


QED Spectral Tests 


Both Perturbative and Nonperturbative 


Perturbative Only 




Static Limit (m2 — >■ oo) 


Reduces to One-Body Dirac Eq. 


Does Not Reduce to One-Body Dirac Eq. 



Table 21: Comparison of TBDE with EEC 



plausible fashion (see Eq. (12) in fT). In the case of the TBDE, natural smoothing mechanisms appear in the 
Dirac formalism allowing one to avoid these adhoc assumptions [28] , [29] El. In addition, in the approach of the 
TBDE, strictly nonperturbative (i.e. numerical techniques) are used for the spectral evaluations. This does not 
appear to be the case for the work of EFG, particularly for the heavy mesons, where use of perturbation theory 
is required because of singular potentials. It may very well be that their adhoc substitutions used in the later 
paper ^Tj will render the use of perturbation theory unnecessary for those mesons. But in that case clear tests 
must ensure that not only do the adhoc substitutions give the same results as the perturbative treatments, but 
that this hold in the sensitive testing grounds of QED [4] , [29] ground states and those of related field theories 
Both approaches display chiral symmetry breaking through the appearance of quark masses. It has been 
demonstrated in the TBDE, however, that the pion mass vanishes in the limit in which the quark mass vanishes 
[2] . That is not demonstrated in the EFG formalism nor for any other potential model formalism that we know 
of for the mesons (In an exception, Sazdjian has demonstrated this using pseudoscalar interactions for the TBDE 

The final point we want to make about the differences is that the wave equation arising from the TBDE has 
been tested in perturbative QED and related field theories. Todorov, and others [10], [36] and [51] showed using 
perturbative methods how their local version of the quasipotential equation displays the accepted QED spectral 
results through order a'^ for two oppositely charged particles with arbitrary mass ratios. The work by Crater 
and Van Alstine [34] and others [29] go beyond this and show that not only do the TBDE display the correct fine 
and hyperfine spectral results when treated perturbatively, but those same results can be recovered when the 
equations are treated nonperturbatively. In [52], a local quasipotential equation closely related to that used by 
EFG in the meson spectrum was shown to also reproduce perturbatively the spectral results through order a* of 
QED for two oppositely charged particles with arbitrary mass ratios. However, the important nonperturbative 

^■^In 1331 . Crater, Yoon, and Wong described some unusual singularity structures of the effective potentials and wave functions that 
show up in Eq. ilfit for singlet and triplet states, in both QED and QCD. In these cases the TBDE lead to effective potentials and 
wave functions that are nevertheless not singular. The most noteworthy case was for the coupled '^Si-^Z?i triplet system, when the 
tensor coupling is properly taken into account. There it was shown that including the tensor coupling is essential in order that the 
effective potentials and wave functions are well behaved at short distances, with the 5-state and D-state wave functions becoming 
simply proportional to each other at short distance (see Appendix C). 

H. Sazdjian has considered chiral symmetry and its breaking in the context of a closely related version of the TBDE. This was 
later discussed by Crater and Van Alstine, [13] . In particular, it is found that the matrix element of the divergence of the axial 
vector current is proportional to the quark mass. This demsonstrates that the quark masses in the TBDE play the same role as 
the quark masses in QFT. In addition, Sazdjian shows that the pion decay constant in the context of the TBDE does not vanish in 
the limit of m,q — > 0. Sazdjian shows analytically in the context of a pseudoscalar confining potential, the existence of a massless 
pseudoscalar meson when rriq — > 0. These are two of the main effects of the spontaneous breakdown of chiral symmetry. In our 
earlier work with Van Alstine |13j . we showed numerically for the case of scalar confinining intereractions that the calculated pion 
mass tends to zero as mq — > 0. In a later work ([2]) we showed, however, that the behavior is not of the square-root relation 
{rrin ~ ^mq/ F-jr). The same behavior appears to hold with the present calculations. 
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tests as done in [55] of the bound state formalism for QED has not been carried out with the local quasipotential 
equation of EFG. 

It is our contention that any relativistic potential model that includes four-vector interactions should, when 
the vector interaction is replaced by its QED counterpart, and confining potentials are set — 0, reproduce the 
the standard QED spectral results. There are two models for which we carry out this test. We limit our test 
to the singlet positronium ground state. The fundamental question we ask is do these two approaches, [Bl and 
[7] , which are used quite successfully for meson spectroscopy give the correct spectral results when restricted to 
QED. The first one we examine is that of Ebert, Faustov and Galkin [7]. 



6.2 Positronium Ground State Spectral Test of the Quasipotential Equation of 
Ebert, Faustov and Galkin. 

The effective Schrodinger equation of this approach restricted to an equal mass bound system for vector in- 
teractions is given in Eqs. (1-4), (13-22) of 7 . For that restriction we have, in the notation of the present 
paper, 





r) ^ Wj 




ei£2 w 
w 4' 


b\w) 


= -(w^-4m^), 
4 ^ ' 


V{r) 


, , N r ,W — \ 

= A{t)[1 + { 

w 


A 


a 



w{w + 2m) iw^ 



(48) 



The unperturbed effective Hamiltonian is 



and the perturbation is for singlet states 

= ^C^-^ f- - - ^)^^(r). (50) 

w r w(w + Zm) w 

Comparing with the nonrelativistic hydrogenic Schrodinger equation 



with ground state energy 



|^-7)V' = fV'., (51) 



(52) 



we see that for the eigenvalue equation (1481) . the total cm. invariant energy from Hq is determined by analogy 
with Eqs. ^ [52]) from 



Thus with 



(53) 



- 4m2 = -A^i^a^ = ^-^'^^ (5^) 
we find that 

^^2m-^ + -^. (55) 
Substituting this into Eq. ([50]) at the appropriate order gives 



H, ^ _(:z:^)-- + :::i^^(r) 

8to r m'' 

ic3 



,ma 2^ 



^S-%r), (56) 
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with their spin-spin terms partiahy cancehng their spin independent contact (Darwin) term while the first term 
is of higher order. The ground state unperturbed wave function is 

exp(-r/ae//) 



(™e//) 



1/2 



Hfi^a ma 

The expectation value is 



1 2 



4 

ma 



Hi) = ——rT3 = (58) 



2 /- 2 



— ) 

mo: / 



and so Eq. (|53p . including (i^i) then becomes, to the appropriate order 







ma'^ 




2 


' 8 


— 4m^ = 




2 



(59) 



so that 



ma 



^ Sma^ ma* ma^ lima* 



This result is in disagreement with the accepted fine structure result of 

ma^ 21ma* , 

w ^ 2m . (61) 

4 64 ^ ^ 

In the Two-Body Dirac equation this spectrum results from an exact solution of the Schrodinger-like form i34j . 

PTHl 17m [24| 



{p^ + 2e,,A~ A^)il; = b'^'4>, (62) 



which yields 



2 + 2/ 



ma^ 21ma* , , 

'"^-^ 6^ + -' (^^^ 



obtained by steps similar to those outlined in Ea.(|5(I)). 

The local quasipotential approach of |7| does not include the term ~ —o? jr^ which gauge invariance 
considerations would demand. It is of interest that if their approach includes this terms then the added 
contribution from this potential is 

^ I ^ 1 , Q?^, ma* 

Y,^-l^^-m^-^^--^- ^''^ 
This, together with the fact that combining this with the earlier results gives the correct added 0(a*) correction, 

lima* ma* 21ma* 
-64 2- = 6^' 

points strongly to the lack of this term as being the cause of the incorrect QED spectral prediction in this 
approach. 

We should emphasize that the closely related formalism of [52] does produce the correct — relativistic 
correction. The difference between the formalism of EFG and that of [55] is that the former does not include 
two loop and iterated Born diagram^ contained in the latter. Those combined diagrams to lowest order in a do 
produce the —A^ = — a^/r^ term which would account for the spectral difference seen in Eq. (|65p (see also |51j). 
Since that —A^ term is not included in the EFG meson spectral formalism it is likely that possibly important 
relativistic corrections for their meson spectrum will be missing. In a private communication, Faustov stated 
that the reason they did not include the contributions of two and more gluon exchange diagrams within QCD 
in calculations of the meson spectra, is that the effects of these diagrams would be contained in the confining , 
long range potential, the origin of which is not known and which is thus added phenomenologically. However the 
—A^ contribution from those two-gluon exchange diagrams due to the Coulomb-like potential A is short range 
and therefore would not by itself contribute to the confining potential. In other words we claim that since its 
effects are short range, it should be considered apart from the phenomenologically added confining interaction. 



^*This equation and its gauge structure can also be seen to result from the equal mass singlet equation version of I I35I1 for S = V = 0, 
A = —a/r or its radial version given by Eq. IIB.Sll l. (It is noted that under these conditions, = 3"I>ss,and ^SOD = ^SOX = 0) 

^^The authors are grateful to Professor R. N. Faustov for pointing out to us the results of 1521 and for their reason that the bound 
state equation used in (7) did not include the —A^ term. 
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6.3 Positronium Test of the Approach of Godfrey and Isgur. 

Although Crater and Van Alstine carried out an earher comparison [2] with this approach , in hght of the 
problem with the above quasipotential approach it is instructive to include a parallel perturbative treatment on 
their different quasipotential equation. Their relativistic Schrodinger equation (see their Eqs. (1-4)) relevant 
for the case considered here has the Hamiltonian which include the spin-spin term in addition to the modified 
Coulomb term (see their Eq. (A15) El). Their equation was of the quasipotential type given in Eq. (IT7|) 
extended to include spin. Here we consider its semirelativistic expansion. 

a 2 „ „ o . a ^ n 1 - 



H - 2vV+^-- + :r^Si-S2V2A--^{p2,-} 
r 6m'^ Ira r 

2m + Ha + Hi, 



Ho 



9 



m r 



(p2)- , 27ra 3 a . ^ 1- 



The ground state unperturbed wave function is 

m = 



exp(— r/a) 

(ttc^) 
2 



1/2 



ma 



(67) 



We find that 



4™ m m m^ 2m^ r 



and with £— ^ms— have 



- {^!{E+-f^) + —-~—^{^{Y>\-m, 69 
4m r 4 2m^ r 



1 



— 



4™ 64 
it , a , a , , ma 

2m^ r ' 8 ^ ' 16 ' 



- — = -1 34^/ c;^exp(-2r/a) = (70) 

" 4m Tra-^ 







and using 



2to^ r 2m2 r r 

£ , ^ a ^ , I , ^ ^ , ma^ ma^ 



-(*-*)--(* — *) = ^ (71) 

m r m 8 2 



and find that 



and so 



ma ,1111, , , 

ma^ ISma^ 

u; = 2m+(/7o) + (i^i) = 2m-— -— . 73 

4 64 

This also does not agree with the accepted result of Eq. (I5T|) . Since the addition of the —a^jr^ term would 
drive this to the other side of the accepted value it is not clear where the error is in this approach. 



^^In their appendix A Godfrey and Isgur modify the Coulomb and contact spin-spin term used here with smearing functions and 
extra non-local parts in order to account for the off mass shell effects not present in the on shell scattering amplitudes from which 
they extract their potentials. We do not include the effects of the gaussian smoothing factors in the determination of the modification 
of the Coulomb term from their Eq. (A15). 
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6.4 Comparison of Two Approaches to Dirac Equation in the Static Limit 



The dynamics of the heavy-hght qQ bound states, particularly the ub and db, should be well approximated by 
the ordinary one-body Dirac equation. In the limit when say m2 — oo details outlined in Appendix A. 5 show 
that our TBDE reduce to the single particle Dirac equation for a spin-one-half particle in an external scalar and 
vector potential, 

(7 • p-/3(e - A) + TO + 5)^- = 0, (74) 

in which e is the total energy of the single particle of mass to. In this same limit Eq. ([35]) (see Appendix A5 ) 
becomes 

{^^ + 2mS+S^ + 2eA-A^ + - 7 + 7 7^ 7 + ^ 7 

2m+S+e-A A\m+S+e-Aj m+S+e-A r 

- (e'-m^)^^, (75) 

which agrees with the Pauli reduction of the Dirac equation (I74p for a single particle in an external scalar and 
vector potential when the first order momentum terms are scaled away (see for example [33] ) . The wave function 
Tpj^ is the upper two-component spinor. From the point of view of the single particle Dirac equation the quadratic 
5'^ and —A^ terms above are not put in by hand but arise naturally from the Pauli reduction. 

We compare Eq. (175]) with the corresponding equations from the two quasipotential approaches, including 
scalar and vector intcractionj^. Referring to Eqs. (1-4), (13-22) of [7] we have the TO2 — >■ 00 limit of 

(p2 + 2e{A + S) + \^ + ^--^^)^ ^ [e^ - to^)^, (76) 
2 77i + £ TO + e r 

in which we have used [i^ — > e. It is evident that for the vector interaction alone (5 = 0), this equation will 
not yield a spectrum perturbatively equivalent to the Dirac spectrum for A = —a/r and for the same reason 
as with positronium, that is, the lacking of the ~A^ term. This would result in incorrect fine structure for 
hydrogen-like atoms (see discussions and footnotes below Eq. (|65p for the reason for this omission). Also, there 
are three parts for the scalar potential interaction that differ from the Dirac equation: the appearance of 2eS 
instead of the expected 2toS', the lacking of the term, and the absence of any scalar Darwin term. The 
appearance of e instead of to can be traced to the use of a common reduced mass for multiplying both vector 
and scalar interactions. Beyond that is the absence of the potential energy terms in the denominators of the 
spin-orbit and Darwin terms. The authors correct this by hand, but their correction does not match the forms 
in the Dirac equation. On the other hand those potential energy terms in the denominators of the Darwin and 
spin-orbit terms of Eq. (1751) provide a natural smoothing mechanism that eliminates such singular potentials 
as delta functions and potentials. For example take the case of S" = 0, A — —a/r. The Laplacian term 
would produce 47ri5'^(r). However, the A term in the denominator would then be evaluated at the origin and 
completely cancel the effects of the delta function. Its perturbative effects are reproduced by the adjacent 3/4 
term. Similarly the behavior of the spin-orbit weak potential form in which the A in the denominator is 
ignored is modified to very near the origin to a less singular 1 /r^ potential by the effect of the A term in the 
denominator as well as the 3/4 term. Similar smoothing mechanics naturally built in to the Pauli structure of 
the Dirac equation occur in the Pauli reduction of the TBDE ( see [28] , [32] , [29] , [33]). 
Referring to A-15, 16 of [6] we have the TO2 — c» limit of 

(v/p^ +m^ + A + S+ '^' ~f ^''^^ )yj = w^j. (77) 

4to^ r 

This Hamiltonian form is missing Darwin terms for not only scalar interactions, but also for vector interactions 
as well. Those are as important for spectral studies as the spin-orbit terms and their lack is a serious defect 
in both these equations. The lack of the vector Darwin term would result in incorrect fine structure for L — 
hydrogen levels. 



7 Conclusion and Future Directions 

The application of Dirac's constraint dynamics applied to the relativistic two-body problem leads quite naturally 
to the Two Body Dirac equations of constraint dynamics when both particles have spin-one-half. This paper 

^^In I53| the static limit Dirac equation was recovered from a two-body quasipotential equation by techniques with some similarity 
to the Gross equation [33]. In that equation the relative energy is constrained by restricting one of the spin-one half particles to 
its positive energy mass shell. This differs from the TBDE which treats both particles off shell but yet constraining the relative 
energy covariantly through P - pip = 0. The EFG equation also has this constraint, but unlike the equation derived in 53., the Gross 
equation and the TBDE, it does not have the Dirac equation as a static limit as seen by a comparison of Eq. 17611 with Eq. 17511 . 
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follows many earlier ones analyzing the structures and applications of those equations. It has several sets of aims 
and results. First we showed that when the interaction structure used in these equations is extended from two 
invariant functions (generated by what we have called A{r) and >S'(r)) to three (not only the two that generate 
an electromagnetic-like interaction and a confining world scalar interaction but also the V{r) that generates a 
confining timelike vector interaction), that the fit to the meson spectrum is improved substantially. However, 
there is still a considerable amount of improvement that is needed, primarily in the radial and orbital excitations 
of the singlet and triplet ground states. Work in progress seeks to extend invariant functions to ones that 
generate covariant pseudoscalar, pseudovector, and tensor interactions. Second, this paper also included 19 
mesons not included in earlier work [2j where only two invariant functions were used. Among those 19 were 
the isoscalar 77 and 77' mesons. Here we developed an approach motivated by some work of Brayshaw [S] which 
introduces a constant symmetric but nondiagonal mass matrix that couples isoscalar qq channels. The three 
parameters introduced are adjusted (when possible) to fit the tt*^, rj and 77' meson masses and then used to 
predict accurately, at least in Th2, the SU{3) pseudoscalar mixing angle . Missing is an attempt to connect this 
mass matrix to the compatibility condition between the constraints Si and S2- 

A glance at Ref. [11] shows that there are no shortages of attempts to stake claims of which relativistic 
two-body truncation of the Bethe-Salpeter equation is most successful. One of the purposes of [5] was to clarify 
some guidelines and important benchmarks that such equations should have when applied to any relativistic two- 
body problem be it for QED bound states, QCD bound states, or nucleon-nucleon scattering. This clarification 
continues in this paper with a fairly detailed analysis of the local quasipotential approach of Ebert, Faustov, and 
Galikin. This approach was chosen among numerous others for two reasons: a) the close connection between 
the minimal dynamical structure of the constraint approach and the early work done on the local quasipotential 
approach by Todorov and his coauthors and b) the extensive phenomenological studies by the local quasipotential 
approach of EFG in meson spectral and decay studies. 

The two quasipotential methods that we discussed have three weaknesses. Neither method uses a wave 
equation to uncover their respective spin-dependent corrections. Rather, they use an on shell version of the 
scattering amplitude and the quasipotential equation for the potential. However, as discussed in both papers, 
they each must make assumptions that allow them to include expected off shell effects. In contrast, the TBDE 
of constraint dynamics include automatically by their mathematically consistent construction, off mass shell 
effects. The second weakness, is that both approaches, when applied to QED, do not produce the correct 
hyperfine structure. This, in our opinion, is a serious but easily correctable drawback in both approaches. 
Let us be precise here about our concerns. In both of these two quasipotential approaches to the relativistic 
QCD potential model, if one turns off the confining interaction and replaces the nonconfining vector potential 
by the Coulomb potential their resultant QED spectrum will be incorrect whether computed perturbatively or 
numerically. This calls into question their QCD spectral results since the resultant relativistic corrections which 
the omitted term —A^ would have contributed is of the same order as the spin-dependent corrections (dependent 
on A' and V^A) which, of course, are not omitted. Even if that correction is made, the wave equation should be 
shown to have (just as can be shown with the one-body Dirac equation) the same spectral results whether treated 
perturbatively or nonperturbatively. In addition, their modeling of the scalar interactions does not conform with 
the approach of the classical and quantum field theories used by Crater, Van Alstine, Yang, Sazdjian and Jolluli 
which supports our choice of vector interaction structurej^. The third weakness is the lack of agreement in the 
m2 — >■ 00 limit with the Pauli-form of the one-body Dirac equation. In particular, the lack of scalar Darwin 
terms in both approaches and vector Darwin terms in the approach of [B] is a serious weakness. 

The constraint approach has been tested against both classical and quantum field theories for both scalar 
and vector interactions. In our construction of the vector potential [32], [24], [29] three primary guideposts were 
used beyond that of a minimal structure. The first is the use of just the barest (lowest order) input from field 
theory, the nonrelativistic Coulomb potential in QED. 

The second is the gauge-like minimal coupling structure Eq. (f22]) of the potential which Todorov postulated 
and was later confirmed in three independent ways: a) Rizov, Todorov, and Aneva |51| demonstrated how the 
gauge structure (particularly the form (e^ — A)'^ — ) arises in perturbation theory at higher order than the 
Born approximation F^. b) by a comparison of the Fokker- Tetrode classical field theory 0(l/c^) expansion for 
the Hamiltonian with the general quasipotential structure of + ^(r,^) — (e^ — m^) [25^ Ej, and c) JoUouli 
and Sazdjian |26j who found a similar structure from nonpcrturbativc quantum field theoretic arguments both 
for scalar and vector interactions. These three arguments demonstrate that the Born approximation structures 

^*In particular, their approach does not include the quadratic structure 2mwS + 5^ implied by the those authors' approaches to 
scalar field theories. 

^^In particular, the quasipotential equation Eq. JTSj for the potential to be used in Eq. II19I I must be iterated to second order 
= TiGTi - T2). It is remarkable that this this gauge like structure postulate Eq. I I22I I anticipates the systematic inclusion of 
higher order terms by the quasipotential formalism. 

^"In this, not only does the gauge like minimal structure (e„ — A)^ appear to be a natural outgrowth of classical 0(l/c^) expansion 
to at least order but also a minimal scalar interaction structure appears so that combined they yield {cw — A)"^ — (m^ -|- S)"^ , 

again, at least through order l/c*. 
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(particularly the first term of (PE)) used to model the QCD potentials in [B]) cannot possibly yield that gauge 
structure and thus cannot yield the correct positronium spectrum when applied to QED. The higher order 
structures in [55] and [5T] as well as the nonlinear hyperbohc structures in Eqs. (|A.2[ IA.3T IA.71 rAr22IIA.25l) argue 
that the Born structure of the first term in (|26p are insufficient and must be supplemented by other invariant 
couplings. The papers by Sazdjian and collaborators [3S], [ID], [55] demonstrate this explicitly, showing that 
pseudovector coupling is essential if done perturbatively in order to get the Dirac equation into an external field 
form in which the minimal structure can be demonstrated. The work in [37], [38) as well as by Sazdjian and 
collaborators show that the vector coupling (see Eq. (|A.10I) ) alone will, when placed in the nonlinear context 
of the hyperbolic parametrization given in Eqs. (jA.2IA.3|) . yield that external field form in which the minimal 
structure can be demonstrated. 

The third guidepost is the use of the relativistic reduced nicLSS TTinj a-nd energy (see Eqs. © and ©). In 
addition to the discussions given in |27| . their appearance in the forms 2mwS + 2ewA as part of the minimal 
structure is also dictated by the same field theory mechanisms discussed in item b) and c) above. Note that 
even though in Eq. (|19|) there is the appearance of another relativistic reduced mass (eie2/tf), in working out 
the quasipotential using the spinors as in Eq. ()25|) that reduced mass is replaced by ruw in the scalar and Sw in 
the vector case. This replacement does not appear in the work of [7]. 

So, while the works of Ebert, Faustov, and Galkin, and Godfrey and Isgur are quite impressive in terms of 
spectral and decay agreements, it would be of value for adherents of those approaches to consider the criticisms 
presented in this paper, related to the ability of those equations to reproduce the static limit Dirac equations 
structures (which would include the —A^ term) and general short range structure related to the —A^ gauge term. 
These criticisms are not about their choice of QCD inspired potentials, but rather about how their relativistic 
wave equations translate the physics of those potentials into spectral results. These thus concern the impact on 
their QCD meson spectral results those two approaches would have based on their field theory connections to 
QED bound states . 

Finally a few words about future direction lines of research related to this paper. The approach given in our 
paper views the meson as a two-body bound state in a first quantized formalism. In place of the nonrelativistic 
Schrodinger equation for two interacting particles, we use the TBDE of constraint dynamics. There are system- 
atic corrections that should follow the completion of this first step which would end with the inclusion of covariant 
pseudoscalar, pseudovector, and tensor interactions, in addition to the scalar and vector interactions we have 
included in this paper. First is a second quantized version the TBDE similar that what has been accomplished 
for nonrelativistic second quantized formalisms by the Cornell group [55] . Tornqvist [56] . and more recently by 
Barnes and Swanson [ST]. The latter includes a microscopic theory [35] of the ^Pq model which includes pair 
production. The aim, as in their recent paper, would be to gain a measure of the effects of two-body meson 
decays on the observed rest mass of the decaying meson. Beyond that would be many body formalisms which 
view the meson more generally as a linear combination of qq -\- qqg -\- ... 

The primary weakness in our results of this paper are: 1) the radial and orbital excitations of the old meson 
spectroscopy; 2) the less than good hyperfine splittings of the cc and bb families compared with the good splitting 
results we obtained with all the other families, including the lightest and most highly relativistic ud ; 3) the 
failure to account for the light {^Pq) scalar mesons; 4) Failure to reproduce the square root Gell-Mann-Oakes- 
Renner relation. It is too early to say if the completion of the first quantization program will rectify any of these 
problems. We point out, however, that should the weakness of our model, relating to the radial excitations of the 
light qq bound states be substantially improved by including pseudoscalar, pseudovector, and tensor interactions, 
this would offer the opportunity of allowing the hypothesis discussed in footnote 13 to be actively considered. 
Whether including these other interactions lead to substantial improvements or not, it will be essential to follow 
in parallel with our relativistic formalism, the second quantized nonrelativistic formalisms developed earlier. 



A Relativistic Schrodinger Equation Details 

A.l Connections Between the TBDE and Eq. ([35]) and Forms for Af, Si in Terms 

of the Invariants A{r),V{r), and S{r) 

Here we present an outline of some details of Eq. ([T4| and its Pauli-Schrodinger reduction given in full elsewhere 
(see [24l [37l 138] |42] ) . Each of the two Dirac equations in (fT4]) has a form similar to a single particle Dirac equation 
in an external four-vector and scalar potential but here acting on a 16 component wave function 5" which is the 
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product of an external part (being a plane wave eigenstate of P) multiplying the internal wave function ip 



4. 

The four ip^ are each four-component spinor wave functions. To obtain the actual general spin-dependent 
forms of those A'^jSi potentials which were required by the compatibility condition [Si,S2]ijj — was a most 
perplexing problem, involving the discovery of underlying supersymmetries in the case of scalar and timelike 
vector interactions |28) , |24| . Extending those external potential forms to more general covariant interactions 
necessitated an entirely different approach leading to what is called the hyperbolic form of the TBDE. The most 
general hyperbolic form for compatible TBDE is 

Slip = (cosh(A)Si + sinh(A)S2)7/' = 0, 

S2ip = (cosh(A)S2 +sinh(A)Si)i/. = 0, (A.2) 

where A represents any invariant interaction singly or in combination. It has a matrix structure in addition to 
coordinate dependence. Depending on that matrix structure we have either vector, scalar or more general tensor 
interactions [37) . The operators Si and S2 are auxiliary constraints satisfying 

Slip = [Sio cosh(A) -t- 520 sinh(A) )V = 0, 

S21P = (520 cosh(A) + Sio sinh(A) )?/' = 0, (A.3) 
in which the Sio are the free Dirac operators 

1 

This, in turn leads to the two compatibility conditions [ST] [37l [39] 

[5i,52]V' = 0, (A.5) 

and 

[Si,S2]V^ = 0, (A.6) 

provided that A(a;) = A(a;j_). These compatibility conditions do not restrict the gamma matrix structure of 
A. That matrix structure is determined by the type of vertex- vertex structure we wish to incorporate in the 
interaction. The three types of invariant interactions A that we use in this paper are 

^c{x±) = -I1I2 — ^^1' ^1 = ~75i752> scalar, 

^j{x±) = /^i/^j^-^^^Ci, timelike vector, 
Q (x ) 

^gix±) — 7i± • 72_L — 2 — spacelike vector, (A. 7) 

where 

1 2 

'~Vr- ■ T ■ T ■ '~V ■ 'y- 

/3, = -7. -A ^ = 1,2. (A.8) 
For For general independent scalar, timelike vector, and spacelike vector interactions we have 

A(xj.) = A£ + Aj + Ag. (A.9) 

The special case of an electromagnetic-like interaction (in the Feynman gauge) corresponds to = or 

Aj + Ag = A£^ = (-7i.P72.P-|-7i^-72^)^^Oi 

= li-l^^^Oi. (A.IO) 
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Our Thl corresponds to a scalar and electromagnetic interaction, 

A{xi_) ^ Ac + Asm- (A.ll) 
Our Tli2 corresponds to a modification of the timelike portion of Asm to 

Aix^) = Ac+Aj + Ag = {-lil2C{x^)+l3^P^Jix^)+-f^^--f2^g{xi_))^, 

-Q ■ (A.12) 

This leads tcl!!! [371 |3H] = d/dxt") 

5iV = ( - G/3iSi • V2 + i?i/3i75i + ^^1751 - G^S2 • d{Cp^-jp,h,,jr^^)i^ = 0, 

= (G/32S2 • Vi + + M2752 + G^Si • a(/:/3i-^/32)75i752)^ = 0, (A.13) 

with 

G ~ expQ, 

V, = Pi -^E,- 9^1],. (A.14) 

The connections between what we call the vertex invariants £, J ^ Q and the mass and energy potentials M;, Ei 
are found to be 

Ml = mi cosh£ +TO2sinh£, 
M2 — m2 cosh£ + mi sinh>C, 
El ~ Si cosh J' +£2sinhJ', 

E2 £2 cosh J' + £1 sinh J'. (A. 15) 

Eq. (jA.13p depends on the standard Pauli-Dirac representation of gamma matrices in block forms (see Eq. (2.28) 
in [5] for their explicit forms) and whcrj^ 

= 75./3z7±.- (A.16) 
Comparing Eq. (IA.13|) with Eq. (|14l) we find that the spin-dependent vector interactions of Eq. (ITil) are [53] 

i^" = ((£1 - El) - z| (72 • dj) 72 • P)P^ + (1 - G)p^ - '-dG ■ 727^^, 

A^2 = ((£2 - E2) + (7i • dj))^i ■ P)P^ - (1 - G)p^l + '-dG ■ 717^^. (A.17) 

Note that the first portion of the vector potentials is timelike (parallel to P'') while the next portion is spacelike 
(perpendicular to P^). The spin-dependent scalar potentials Si are 

51 - Ml - mi - ^G72 • dC, 

52 = M2 - m2 + ^G7i • dC. (A.18) 



Eq. (|A.17P simplifies to 



At = ((£1 - Ei))P^ + (1 - G)p'i - ^aG • 727^, 
A'i = ((£2 - i?2))P^ - (1 - G)p^^ + U)G ■ 7i7^ (A.19) 
for electromagnetic-like interactions. 

^^In short, one inserts Eq. l lA.31 1 into I IA.211 and brings the free Dirac operator HA. 411 to the right of the matrix hyperboHc functions. 
Using commutators and cosh^ A — sinh^ A = 1 one arrives at Eq. IIA.131 1. The structure of these equations are very much the same 
as that of a Dirac equation for each of the two particles, with Mi and Ei playing the roles that m + S and e — A do in the single 
particle Dirac equation {ZH- Over and above the usual kinetic part, the spin-dependent modifications involving CPi and the last 
set of derivative terms are two-body recoil effects essential for the compatibility (consistency) of the two equations 

Just as x^^ is a four vector, so are 7^ (in the sense of Dirac) and Pf. Thus, the matrix structures of the time-like and space-like 
interactions in Eq. IIA.71 1 are 7^72 a^nd Ti ' T2 oi^ty ir' the cm. system due to the fact that from Eq. IIA.8I I. /3j = 7^ only in the cm. 
frame. Likewise, = (0, S) only in the cm. frame just as is x'l = (0, r) in that frame only. 
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We have chosen a parametrization for £, J , and Q that takes advantage of the Todorov effective external 
potential forms and at the same time will display the correct static limit form for the Pauli reduction (see Eq. 
([75]) V The choice for these parametrizations is fixeded due to the fact that for classical [35] or quantum field 
theories |26] for separate scalar and vector interactions the spin independent part of the quasipotential $^ 
involves the difference of squares of the invariant mass and energy potentials (M^ and Ei respectively) 

Ml = TO^ + Im^S + 5*2; El = e\ - 2e^V + , (A.20) 

so that 

Mf - E^ ^ 2m^S + S^ + 2e.^V - - b^{w). (A.21) 

Strictly speaking, the forms in Eq. (jA.20[) and Eq. (|A.21I) are for scalar and timelike vector interactions. 
Eqs. (ITi|) and (IA.13|) involve combined scalar, electromagnetic-like, and separate timelike vector interactions. 
Without the separate timelike interactions this amounts to working in the Feynman gauge with the simplest 
relation between space- and timelike parts, (see Eqs. (jA.lOl) . (jA.lip . and [13l[2]). In the general case the mass 
and energy potentials in place of Eq. (|A.20p are respectively 

Mf = mf + exp{2g){2m^S+S^), (A.22) 
E^ = exp{2g{A))i{e^~Af-2ey,V + V^), (A.23) 

so that from Eq. ()A.15|) . 



.,, ^ml + exp(2g)(2m^5+52) + y/mj + expi2g)i2m^S+S^) 

exp(£) = exjp[C{S,A))^ ■ , (A.24) 

mi -I- TO2 

exp(J) = exp(J(F,yl)) = exp(^)-!^ , 

El + £2 

with 

exp(2g(^)) ^ (^3^17^ - (A-25) 
Below we present, the consequent connections to the invariant interaction functions A, l/,and S. 

a) In the case of electromagnetic interactions {V — 0) with scalar confinement (Thl), we have 

1 E-\ + Eo 

J = -g= log{l-2A/w)^\og^—^, 
I ui 

Ef = exp{2g){e^~Ay, 

= mf + ejip{2g) {2m^S + S^) , (A.26) 

and the spin-independent (SI) minimal coupling comes from 

exp(2g)p2 + - £f ^ p2^exp(-2g)(Mf -£;2) 

= + ^si - (A.27) 

and appears as 

- 6^ = 2m^S + 52 + TO^(1 - 2A/w) - (e, - A^ = 2m^S + + 2e^A -~ A^ - \? . (A.28) 

b) In the case of pure timelike vector interactions and scalar interactions (with no electromagnetic-like inter- 
actions this model is not appropriate for meson spectroscopy) we have 

rr rr 1 ^10 + -^20 

J Jo = log , 

w 

El ^ ei-2£^V + V\ 
a = 0, (A.29) 

and the spin- independent minimal coupling appears as 

^si = 'im^S + S'^ + 2e^V -V"^. (A.30) 
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c) When we include independent timelike and electromagnetic-like simultaneously together with scalar inter- 
actions (Th2) then we have 

-g = ^log{l-2A/w), 

J = log , 

w 

Ef = exp(2^)( (e, - A)^ ~ 2e^V + V^) = exp{2g){El - 2e,A + A'') (A.31) 
and the spin-independent minimal coupling appears like 

$s/ = 2m„5 + S^ + 2£^A - + 2s^V - V"^. (A.32) 

A. 2 Details on Eq. (I35l) 

The Klein-Gordon like potential energy terms appearing at the beginning of the Pauli form (j35]) arise from 

m2 _ Ef = exp(2^;)[2TO„S' + + 2eyjA - A^ + 2£^V -V'^ - \?{w)\. 

To obtain the symbolic Pauli form of Eq. and the subsequent detailed form in Eq. (1551) involves steps 
similar to those used in the Pauli reduction of the single particle Dirac equation [33 but with the combinations 
— ± ip^ and x± = V'2 i "03 instead of the individual single particle wave function. This reduces the Pauli 
forms to 4 uncoupled 4 component relativistic Schrodinger equations [iQl [HI |38l [2] |42]. We work in the cm. 
frame in which P = (1,0) and f — (0,f). The final four-component wave functions 'ipj.,rij_ that appear in Eq. 
([55]) are defined by |^ 

(j)^ = exp(J" + /C(Ti-f(T2-f)?/'^ = (exp J')(cosh/C + sinh/Ccri-fcr2-f)7/;_|_, 

X± = exp(J^ + /C(Ti-f(T2-f)?7_|_ — (exp J^)(cosh/C -I- sinh/C(Ti-f(T2-f)77_|_, (A. 33) 

in which 

1 V 

T = -iog g, 

2 £2^1 + ei"i2 
V^E2Mi + E1M2, 

^^(W)^ (A.34) 

In analogy to what occurs in the decoupled form of the Schrodinger equation for the individual single particle 
wave function, this substitution has the convenient property that in the resultant bound state equation, the 
coefficients of the first order relative momentum terms vanish. 

Using the results in and [33] we obtain for the general case of unequal masses the relativistic Schrodinger 
equation (|35p that is a detailed cm. form of Eq. (I16p . In that equation we have introduced the abbreviations 



2(J-' + l/r)(cosh2/C-l) ^,2 2/C'sinh2/C 

$D = ^ + ^ +)C H V^F + m{r), 

r r 

J" (J'' + l/r)(cosh2/C-l) /C'sinh2/C 

'^so — 1 , 

r r r 

*soD ^{l' cosh 2/C - q' sinh 2/C) , 
= (q' cosh 2K.-1' sinh 2/C) , 

2/C'sinh2/C 2(J-'-f 1/r) (cosh 2/C- 1) 2T'IC' V^JC 



^ss = k{r) 



3r 3r 3 3' 



where 



1^ , , , (3J-'-/C^ + 3/r) sinh 2/C ^ (J"' - 3/C^ + l/r)(cosh 2/C - 1) , . 

$T = T:\n[r) H \ V2T K, - V /C , 

6 r r 

^,cosh2/C-l /C' (J"' + l/r)sinh2/C 

$soT = -/C' + ^ — , (A.35) 

r r r 

fc(»^) = 3V (/c + e^) ^ -^^ , 

n{r) = 1[V^/C - iv^^? -f + ^'(g' - 2/C')], 

3 2 2r 

m(r) = -iv2g+ + g'^' - /C'2, (A.36) 
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and 



1 E2M2 - EiMi , , 

I'ir) - p\r\ p\r + ^ ' 
2r tj2Mi + tjiM2 

1 E2Mi-E^ 

(The prime symbol stands for d/dr, and the exphcit forms of the derivatives are given in Eq. (jA.38p befow). For 
L — J states, the hyperbohc terms cancel and the spin-orbit difference terms in general produce spin mixing 
except for equal masses or J = 0. For ease of use we have listed in Appendix A3 the explicit forms that appear 
in the above $s in Eqs. (|A.35p - (|A.36p in terms of the general invariant potentials A{r)^ and S{r). The 

radial components of Eq. p5l) are given in Appendix B. 

A. 3 Explicit Expressions for terms in the Relativistic Schrodinger Equation ( I35p 
from A{r),V{r) and S{r) 



Given the functions Air) , Vir), and S{r) for the interaction, users of the relativistic Schrodinger equation ([55]) 
will find it convenient to have an explicit expression in an order that would be useful for programing the terms in 
the associated equation (|A.35P - (|A.37p . We use the definitions above given in Eqs. ()A.22h - (|A.25p . and ()A.34p . In 
order that the terms in Eq. (jA.35P be reduced to expressions involving just A{r), V{r), S{r) and their derivatives, 
we list the following formulae 

, ^ {C + J'){E2M2 + EiMi) 

2{E2Mi + E1M2) 
A' 

c = 

W-2A' 

^ _ M[ _M'2 _ w fS'{m^ + S) , {2m^S + S^)A'^ 
J' 



Ah Ml M1M2 V W-2A (w-2Ay 

- A){e2 -A)+ 2e^V - V^] + (e^ - V)V') 
E2 El EiE2{w-2A)/w 



Also needed are 



X' = iSL_ll_ (A.38) 



cosh2^ = 1 (£1 + £2)(A^i + M2) , (mi + m2)(^i + E2) 



2 \{mi + m2){Ei + E2) (ei + e2)(Mi + M2) 

2 \{mi + m2){Ei + E2) (ei + e2)(Mi + M2) ' ' ^ ' 



and 



^ {V^C + V^J){E2M2+EiAh) ,,,^,.2 {^l-^lf 

2{E2Ah+EiM2) ^ ' 2{E2Mi+EiM2f 



M1M2 

w fV^S{m^ + S) + S'^ 4S'{m^, + S)A' + {2m^S + S^)V^A 4{2m^S + S^)A''- 



M1M2 V W-2A {w- 2AY (w - 2A)3 

V^J = -[^il±Ii)j'-2g']J'-''-^l^{V^g[{ei-A){e2-A) + 2e^V-V^] + {e^^V)V^V 
-g'^iw - 2A)2 - V'^ + 2V'g'{e^ - V)} 

v^g - ^ +2g'\ (A.40) 

w — ZA 

The expressions for fc(r),m(r), and n{r) that appear in Eqs. (jA.35l) ) are given in Eqs. (jA.36p . They can be 
evaluated using the above expressions plus 

V^^ ^ (A.41) 

The only remaining parts of Eq. (IA.35P that need expressing are those for I' and q' given in Eq. (|A.37I) Using 
Eq. (|A.34p they can be obtained in terms of the above formulae. 
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mim2 mi + TO2 
A e^,V A + V 



A. 4 Weak Potential Limits of Quasipotentials 

The weak potential forms of the quasipotentials are needed for working out perturbative spectral corrections. 
For weak potentials, we take Si = rrii wherever they appear in the potentials and assume \V\, \S\ << rrii. 
Thus 

^ m^oS S 

w Si£2 rill + 

A A 

g ^ (A.42) 

W mi + TO2 

Among other results, this will allow us to see how the scalar interaction contributes oppositely to the spin-orbit 
and Darwin terms from both vector interactions to lowest order. In that same limit, the dominant portions of 
$'s in Eq. ^ are 

* V72T- 1 ,V^{S - A-V){ml + ml) ^. 



4r (mi +7712) m2mi 
"){rnl 
\mi 

1 {S' - A' ~V'){ml-ml) 



r 1 AS' - A! -V'){ml + m\) 

^so ^ > 7 A- '-^ - A') 

r 4r (toi + 7712] m^mi 



^SOD ^l' ^ - J , ^ . 

vr \m\ ^ m^) mim2 

<^sox ^ 

3 b[mi+m2) 



T 
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jS + A + V) 
9(mi + m2) ' 



/C' (S' + A' + V) 
r Zr[mi + m2) 

Note how in the Darwin and spin-orbit terms the S' — V' dependence shows how the scalar and timelike confining 
effects tend to cancel for ^ < 1. As anticipated, only the Darwin and spin-orbit terms survive in the static limit 
when one of the two particles becomes very massive. In that limit (say TO2 — ^ 00), the two spin-orbit terms $so 
and $soD combine to 

^ L • (<Ti + CT2) L • (a, - CT2) = -L • cTi . (A.44) 

A. 5 Single Particle Limit (1112 00) of the TBDE 

In addition to using p2 (m-2, 0), Pi = (si, p), the single particle limit of the TBDE is obtained by substituting 
m2 — 00 in the expressions for the various potentials listed in Appendices Al, A2, and A3. To determine that 
limit, we use that the total cm. energy w — £2 + £1 ^ m2 + Si- In that case using the expressions for and 
£w in Eqs. ^ and ([8|) we find 

m-uj —7- mi = m, 

£w £1 = £, 

G = (l-2A/(m2 + e))-i/2^1, (A.45) 
and thus for the quantities related to the mass potentials, we have 

Ml \J{m + Sf = m + S, 



M2 - TO2 Y "^2 + ^wS" + S'2 - m2 0, 

exp £ = exp ■ ) ^- 1, 

mi + m2 

Si S, 



(A.46) 
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For the energy potentials we have, since P (1, 0) that 



El yje^ - 2e{A + V) + A^ + V^ = e-U 



e2-E2 e2-\el-2e2A-2eV + A^ + V^ 



£2 - £2\/l - 2A/e2 A, 

exp(J^) = exp(;^3J-^) ^ 1, 

£1 + £2 

A'i ^ f/(l,0), 

^ ^(1,0). (A.47) 

The TBDE ^ thus reduce to 

(7i-p-/5i(£-C/)+m + 5)V. = 0, 

(-/32(to2 - A) +m2)V' -> -TO2(/32-l)V' = 0, (A.48) 

effectively to a single ordinary one-body Dirac equation (1741) for a particle in combined external scalar and time 
only component vector potential. Note that [/ 7^ A + K Doing the usual Pauli- reduction yields for the upper 
component wave function the Schrodinger-like form which is the same as the m2 — > 00 limit of Eq. (j35]). That 
limit can be readily seen from Eq. (jA.34p and Eqs. (|A.45IIA.47| which shows 

^ 1, m + S + e-U 

T -> - loa 



2 m- 
1 S" - [/' 



2m + S' + £-A' 



1 V^S-V^L/ 1 / S' -U' 



2m + 5' + £-[/ 2\ra^S^e-U 

c',j\g',]c' ^ 0, 



1 

r 



1 i;2M2 - glMi ^ 
^ ^ --2^i?2Mi+i?iM2^^ ^'^^ 

1 S2M2 r 



2r E2M1+ E1M2 r 



(£' + J') ^ . (A.49) 



From the last line we obtain 



2F' S' — U' 
$soL • (cr + (T2) + ^souL • (cr^ - 0-2) = L • cri = — — -L • cri. (A.50) 

r m + 6 + £ — (7 



and hence Eq. (|35|) reduces to Eq. (|75p . In Section 6.4 we display our static limit equations for Thl. 



B Radial Equations 

The following are radial eigenvalue equations, [H], [33], corresponding to Eq. ([35]) . For a general singlet ^ Jj 
wave function wlsj = ^joj = ^^o coupled to a general triplet '^Jj wave function wjij = wi, the wave equation 

{-^ + '^^^^-^ + 2m^5 + 52 ^ 2£,„A -A^ + 2en,V -V^ + $d-3$5s}«o 
+ 2VJ(JTT)($soi5 - $sox)wi 

= fo^^o, (B.51) 

^^Strictly speaking, in the static limit, the square roots forms obtained from Eq. |)A.22|I imply A4i — |m + 5| > 0. For for large 
distances, A — > 0, and so Eq. I IA.2311 implies Bi — > |e — V\ and for short distances V — > 0, i?i — > |e — j4| . These could possibly be in 
opposition to the forms with indefinite sign of m + S and e — or e —A that would appear in the one body Dirac equation. In 
our applications to meson spectroscopy S is always positive so we need not worry about the use of the square root form of Mi in 
the computation of Mi. Since V is positive and increasing with distance, it is possible that at large enough distances e — V < 0. 
In that case use of the positive root for the square root would not agree with the sign of e — y. We found that for the ^Sq bu and 
bd mesons that does occur near the very end of the integration cutoff distance of about 1.92 Fermis. That is so close to the end, 
however, that use of the positive square root is unlikely to have any effect on the spectral results (for the ^P2 bu and bd mesons 
e — V remained positive throughout the integration region.) Future work may address the theoretical problem of how to make the 
approach to the static limit of the square root forms smooth and exact. 
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is coupled to 



- 2$5o + <^ss + 2$T - 2$5ot}wi + 2Vj(J+ l)($soD + *sox)wo 

= 62„i. (B.52) 

For a general 5=1, J = i+1 wave function vj-nj = w+ coupled to a general 5*=!, J = L — 1 wave function 
vj+iij = the equation 

ar r 

+ 2( J - l)$so + *ss + ^7^(*soT - $t)}w+ 



is coupled to 



2J+1 ^ ^^^ ~ ^^'^ + 2)$50T}f- 

(B.53) 



2( J + 2) 

- 2( J + 2)$so + *ss + 2j^/ (^gOT - $t)}w- 



2./J(J+ 1) , 

2J +1 ^ ^^^ + ^^"^ " 1)*sot}i^+ 

(B.54) 



For the uncoupled ^Pq states the single equation is 



2 

- 4$so + ^ss + 4($50T - $t)}w- = (B.55) 

C Two Body Dirac Equations for QED 

The Schrodinger-like form Eq. ((35)) of the TBDE given in Eq. (IT4t can be used for QCD bound state (meson 
spectroscopy) and for QED bound states (positronium, muonium, and hydrogen- like systems). In our meson 
spectroscopy work presented in this paper we use the three invariant functions S(r), A{r), and V{r) specified 
in Sec. 3. Once these are specified then so are the three vertex invariants C{r), G{r), and J'(r). They, in 
turn, fix the quasipotentials given in Eq. (|X35|) and below to (KaH that appear in Eq. p5p and its radial 
forms of Appendix B. To make the transition from QCD meson bound states to QED bound states we simply 
take S{r) — V{r) — and A{r) = ^a/r (we remind the reader that only in the cm. system is the invariant 
r = \f^\_ equal to |r|). Our QED spectral results then follow from solving nonperturbatively (i.e. numerically 
or analytically) the radial eigenvalue equations of Appendix B. For example, since — 3$ss for 5 = ^ = 0, 
the equal mass spin singlet equation (jB.511) collapses to 

^2 J(J + 1) 2e^a ^ 

|-T-2^ 2 T\'"o = bvQ, (C.56) 

ar r r r 

which has the analytic spectral solution given in Eq. for J ~ 0. As long as J(J + 1) — > —1/4, since 
the above equation take the limiting form of 



, J(J+1) a% „ 

{-^+ \,2 -^}^o^O, (C.57) 



it is clear that the effective potential is nonsingular, thus allowing a well defined eigenvalue solutioro. The paper 
[33] demonstrated that the corresponding short distance behavior of the potentials in the other radial equations 



particular, the short distance radial behavior 
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of Appendix B for the other QED bound states also ahow weh defined eigenvalue solutionf[fl. Beyond that, 
numerical solutions of these eigenvalue equations yield spectra (not limited just to the singlet ground states) 
in agreement with standard perturbative 0{a*) results |29] . In electron volts the numerical binding energy for 
the singlet ground state of positronium from Eq. (jB.51|) is -6.8033256279 versus the perturbative result of Eq. 
([gT|) or m(-aV4 - 21a'*/64) = -6.8033256719. The difference in units of is -6.08E-05 which is on the 

order of a^, so that, as expected from (|6ip the difference is on the order of ma^. The corresponding numerical 
binding energy in electron volts for the triplet ground state from the coupled equations Eqs. (|B.53p and (jB.54|) is 
-6.8028426132 versus the perturbative resuU f^S] of m{-a'^/4 + a''/192) = -6.8028426636. These resuhs do not 
include the effects of the annihilation diagram. The difference in units of is 6.97E-05 which is on the order 

of a^, so that the difference is also on the order of ma^. These results (from a very extensive list of numerically 
computed spectral results in [53]) taken together, represent crucial tests of this formalism, ones that have not 
been demonstrated in any other relativistic bound state formalism. In fact, the authors of [4] have found a 
particular quasipotential formalism that does give such agreement, but only for the ground state. They also 
demonstrate that several well know two-body relativistic bound state formalisms (including the Blankenbecler- 
Sugar formalism and the formalism of Gross) fail this test. Let us be explicit about the implications of either the 
failure of this test or the lack of performing this test. When one proposes a new bound state formalism such as 
Dirac did in 1928, it was essential that at the very least it reproduce the nonrelativistic hydrogenic spectrum. 
Beyond that it provides two other remarkable results. First of all, by way of an order 1/c^ expansion, it gave 
an order a'^ perturbative correction to the nonrelativistic spectral results. Remarkably, an exact solution of the 
same equation later by Darwin gave spectral results which when truncated to order agreed precisely with 
the perturbative treatment, and, at the time, with experimental fine structure measurements. Breit, in the 
development of his two-body equation gave us an equation with interactions beyond the Coulomb potential that 
ultimately reproduced, when treated perturbatively, spectral results for two-body systems that agreed through 
terms of order with experiment. Unlike the one-body Dirac equation which has an exact spectral solution 
which agrees, at order a^, with its perturbative solution, the Breit equation has no exact or for that matter 
numerical solution which agrees, at order a*, with its perturbative solution. The same could be said for most all 
other two-body relativistic treatments proposed since then, with the two exceptions (^S| and |3]) noted above. 
One's two-body formalism having an agreement of its nonperturbative treatment with its perturbative treatment 
of the spectra in a well established field theory such as QED, in our opinion, should be regarded as a necessary 
hurdle to pass before going on to apply these formalisms to potential models for meson spectroscopy. 
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